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FOREWORD 


This  report  summarizes  part  of  a  study  performed  by  the  Auburn 
University  Electrical  Engineering  Department  concerning  the  glint 
problem  as  applied  to  radar-guided  homing  missiles. 

This  volume  of  the  final  report  describes  the  unclassified  glint 
literature  and  a  model  used  in  the  study  of  the  problem,  while  Volume 
I  contains  a  discussion  of  prior  classified  glint  investigations  as 
well  as  a  summary  of  other  research  performed  as  a  part  of  the  project 
"Computer  Smoothed  Homing."  This  work  was  performed  under  the  asupices 
of  the  Engineering  Experiment  Station,  Inc,,  Auburn,  Alabama,  and 
was  sponsored  by  the  Future  Missile  Systems  Division,  United  States 
Army,  Department  of  Defense.  Technical  direction  of  this  effort  was 
under  the  Army  Missile  Command,  Redstone  Arsenal,  Huntsville,  Alabama, 
under  contract  DA-01-021-AMC  14525(Z). 
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ABSTRACT 


The  unclassified  glint  literature  is  described  and  discussed, 
with  particular  emphasis  on  applications  of  diversity.  A  computer- 
implemented  model  based  upon  the  concept  of  coherent  summation  of 
contributions  from  independent  scattering  complexes  is  developed.  The 
application  of  the  model  is  illustrated  using  synthesized  scattering 
patterns.  Data  processing  techniques  for  frequency  and  space  diversity 
are  investigated  with  the  model.  The  results  indicate  a  significant 
reduction  in  RMS  glint  error  is  possible  by  the  appropriate  utilization 
of  diversity  techniques. 
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FINAL  REPORT 


CONTRACT  DA-01-G21-AMC  14525 (Z) 

VOLUME  II:  THE  DETERMINATION  OF  GLINT  FOR  A 
COMPLEX  RADAR  TARGET 

R.  J.  Sims  and  E.  R.  Graf 

INTRODUCTION 

As  indicated  in  Volume  I  of  this  final  report,  the  research 
performed  under  this  contract  may  be  divided  roughly  into  three 
phases.  Phase  I  dealt  with  the  design  and  parametric  study  of  a 
dual-mode  sequentially-lobed  antenna  array  and  its  use  in  a  strapped- 
down  radar  seeker  head,  and  also  with  the  radar  signature  of 
vibrating  targets.  Phase  II  included  the  study  of  the  target-image 
effect  and  the  glint  induced  by  simple  geometric  targets  which 
subtend  a  significant  part  of  a  tracking  antenna  beamwidth.  Phase  III 
evolved  from  the  latter  part  of  Phase  II.  The  goal  of  this  effort 
was  to  compile  a  bibliography  of  the  glint  literature,  analyze  these 
prior  investigations  with  particular  emphasis  on  diversity  techniques 
for  reducing  glint  error,  and  to  use  this  analysis  to  guide  the  formula¬ 
tion  of  a  glint  model  for  a  complex  radar  target.  The  results  of 
Phase  I  and  Phase  II  are  summarized  in  Volume  I,  which  also  includes 
a  discussion  of  the  classified  glint  literature  analyzed  under  Phase 
III.  It  is  the  purpose  of  the  present  volume  Co  report  the  unclassified 
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portion  of  the  Phase  III  research  effort. 

Back-scattered  wave-interference  phenomena,  which  occur  whenever 
a  target  of  complex  chape  is  illuminated  by  a  radar  with  radiation 
wavelength  much  less  than  the  physical  dimensions  of  the  target,  can 
give  rise  to  several  different  target  scintillation  effects.  Among 
these  are  fading  of  the  radar  echoes  with  changing  aspect  angle, 
variation  of  the  apparent  range  to  the  target,  variation  of  the 
received  radar  frequency  about  the  Doppler-shifted  value  expected  from 
the  radial  component  of  the  target  velocity  vector,  and  variation  in 
the  apparent  angular  location  of  the  target  with  respect  to  the 
radar.  This  latter  effect  (glint)  can  introduce  errors  of  such  magni¬ 
tude  as  to  indicate  the  angular  location  of  the  target  to  be  many 
target  spans  outside  of  the  region  occupied  by  the  target  itself.  With 
the  increasing  radar  precision  being  demanded  by  the  space  and  missile 
age,  such  errors  may  soon  be  intolerable. 

This  report  develops  a  technique  which,  with  the  aid  of  the 
modern  digital  computer,  may  provide  a  means  of  describing  the  glint 
phenomenon  for  a  complex  radar  target.  In  order  to  demonstrate  the 
basic  glint  behavior,  a  simple  example  will  be  enumerated.  This  will 
be  followed  by  a  review  of  the  techniques  and  results  of  various  glint 
investigations  in  the  form  of  a  brief  history  of  the  study  of  the 
problem.  Finally,  in  the  light  of  the  previous  investigations,  a  glint 
model  will  be  developed.  Emphasis  will  be  placed  upon  maintaining 
sufficient  generality  and  flexibility  in  the  model  to  permit  the  analysis 


3 


of  both  theoretical  target  configurations  and  actual  complex  radar 
targets.  Also}  various  forms  of  parameter  diversity  vi5'i  i;° 


discussed  in  relation  to  their  applicability  in  reducing  glint  errors 
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region  may  be  assumed  to  be  linear  functions  of  0.  (This  is  equivalent 
to  expanding  the  pattern  in  a  Taylor’s  series  about  80  and  retaining 
the  first  two  terms.)  Letting  the  magnitude  of  the  pattern  slope  in 
the  crossover  region  be  p,  the  amplitude  response  may  be  written 


E  =  E0[l  +  p(0  -  eo)] 


(2-1) 


where  F,  Js  the  response  at  boresight.  If  a  single  target  reflector 
is  located  at  an  angle  0.,  one  response  of  lobe  A  is 


E^=  K[1  -  p(0^.  -  0o)  jcos  cat 


(2-2) 


and  of  lobe  B  is 


Eg=  K[1  +  p(3^  -  0Q)]cos  cot 


(2-3) 


where  co  is  the  radian  radar  frequency  and  K  is  a  constant  representing 
such  parameters  as  target  cross  section  and  system  gain.  The  radar 
receiver  derives  an  angular  error  signal  by  sequentially  sampling  and 
square-law  detecting  the  outpucs  of  the  two  antennas,  passing  the 
resulting  voltages  through  a  low-pass  filter  which  rejects  all  but  the 
d-c  and  low-frequency  terns,  and  forming  the  difference 


e  -  eb  -  Ei 


(2-4) 


where  the  bar  indicates  the  filtered  value.  Carrying  out  these  opera¬ 
tions, 


=  K2[l-2p(8T-0o)  +  p2(0T-8o)2]cos2  cot 
=  K2[l-2p(0T-0  )  +  p2(0T-8o)2]tl/2  +  1/2  cos  2u>t] 
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and 

E|=(i/2)K2[i-2p(eT-0o)  +  p2(eT-e0)2]. 

Similarly, 

E2»(1/2)k2[1+2P(0t-8o)  +  p2(0T-eQ)2] 

thus, 

e  =  2K2p(0T-0o)  .  (2-5) 

In  operation,  this  error  signal  is  used  to  command  a  mechanical  posi¬ 
tioner  which  aligns  the  radar  boresight  axis  with  the  target  return. 
Thus,  when  e  is  reduced  to  zero,  0Q  =  0-j.  and  the  radar  is  properly 
tracking  the  target. 

Now  consider  the  case  of  a  target  comprising  two  reflectors, 
located  at  angles  0,^  and  0^2*  For  a  stationary  target,  the  frequency 
of  the  two  resulting  returns  will  be  the  same,  but  the  relative  phase, 
a,  will  depend  upon  range  and  aspect  angle.  Letting  a  be  the  ampli¬ 
tude  ratio  of  the  two  reflectors,  the  lobe  responses  are 

Ea  =  K[l-p(8rj.j_-80) ]cos  wt  +  aK[l-p(0,p2-®o)]cos  (cot+a)  (2-6) 

and 

Eg  =  K[l+p(0Tl-0o)]cos  cot  +  aK[l+p(8T2-0o)  ]cos  (cot+a).  (2-7) 
Again  carrying  out  the  squaring  and  filtering  operations, 


e2  =(i/2)K2[i+2P(eTl-eo)  +  p2(8T1-e0)2] 

+  aK2cosa[l+p(6T1-0o)  +  p(0T2_6o)  +  P2(exl_eo^eT2-6o^ 

+(l/2)a2K2[l+2p(0T2-0o)  +  p2(eT2-60)2]  (2-8) 

where  the  upper  and  lower  signs  correspond  to  lobes  B  and  A  respec¬ 
tively.  The  error  signal  is  then 


e  <=  2K2p(0Tl-0o)  +  2a2K2p(0T2-0o) 

+  2aK2cosa[p(0Tl-8o)  +  p(0T2-0o1 

=  2K2p|8^(l+a  cosa)  +  S^te2**  cosa)-0o(l+a2+2a  cosa)f.(2-9) 


Letting  the  angular  spacing  between  the  reflectors  be  0d3:0t2-0t^  and 
letting  e  =  0  (an  on-target  indication) ,  the  boresight  position  becomes 


[a2  +  a  cos  a 
1  +  a2  +  2a  cos  a 


(2-10) 


An  examination  of  (2-10)  for  various  cases  is  instructive.  If 
reflector  1  is  much  larger  than  reflector  2  (a-K3)  ,  then  Sq-K^I*  an<^> 
conversely,  if  reflector  2  dominates  reflector  1  (a-*00),  ®0‘>®ti+6q:s®x2* 
These  cases  are  equivalent  to  accurate  tracking  of  a  single  reflector 
as  described  above.  If  the  ranges  to  the  reflectors  are  such  that 
the  received  signals  are  in  phase  (<3=0) ,  then  (2-10)  reduces  to 


(2-11) 


and  the  radar  accurately  tracks  the  amplitude  "center  of  gravity"  of 
the  two  reflectors.  However,  some  values  of  a  can  lead  to  erroneous 


tracking.  Consider  the  case  where  the  propagation  paths  from  the 
reflectors  is  such  that  a  =  180°.  Then  (2-10)  becomes 


0o  " 


eT1  +  e 


D  a  -  1 


(2-12) 


If,  for  example,  a  =  0.5,  then  (2-12)  gives  0Q  =  0^-^  -  0^,  and  the 
radar  tracks  a  point  outside  of  the  region  between  the  reflectors. 

Thus  the  boresight  position,  or  balance  point,  is  seen  to  depend  upon 
relative  amplitude,  relative  phase,  and  the  spacing  between  the  two 
reflectors.  This  relationship  is  depicted  in  Fig.  2-3,  which  is  a  plot 
of  (2-10)  with  amplitude  ratio  as  a  parameter.  The  boresight  position, 
0,  is  plotted  in  units  of  the  spacing  between  the  two  reflectors,  , 
when  the  balance  point  is  above  or  below  the  region  defined  by  0=+l/2, 
the  tracking  point  is  outside  of  the  region  between  the  reflectors. 

(Of  course,  physical  limitations,  such  as  antenna  beamwidth,  would 
prohibit  the  very  large  errors  indicated  as  a-KL.) 

While  this  simple  model  should  not  be  considered  as  realistic,  it 
does  demonstrate  how  glint  can  occur.  It  is  also  useful  in  that  it 
suggests  several  possibilities  for  obtaining  more  accurate  tracking 
information,  in  the  presence  of  glint. 

If  bv  some  means  two  or  more  independent  samples  of  the  apparent 
target  position  could  be  obtained  by  the  radar,  then  it  is  reasonable 
to  presume  that  some  computer  processing  or  smoothing  technique  could 
yield  a  probable  value  more  accurate  than  that  indicated  by  a  single 


sample.  Because  the  glint  error  is  strongly  influenced  by  the  relative 
phase  between  the  echoes  from  the  different  parts  of  the  target,  some 


11 


means  of  changing  this  phase  would  permit  independent  sampling  of  several 
values  of  target  position.  The  phase  is  a  function  of  the  difference 
in  propagation  paths  from  the  reflectors,  thus  viewing  the  target  from 
a  different  aspect  angle  would  change  the  relative  phase.  This  could 
be  achieved  by  tracking  (or  illuminating)  the  target  from  several  loca¬ 
tions,  and  might  be  termed  space  diversity.  The  separations  would  be 
chosen  sufficiently  large  to  yield  independent  samples.  The  phase  also 
depends  upon  the  radar  wavelength.  If  the  target  is  illuminated  by 
frequencies  over  a  sufficiently  large  bandwidth,  independent  balance 
points  result.  This  technique  is  called  frequency  diversity.  Finally, 
if  the  phase  of  the  return  from  the  reflectors  is  sensitive  to  the 
polarization  of  the  incident  illumination,  polarization  diversity  might 
yield  the  desired  samples.  It  should  be  noted  that  each  diversity 
technique  can  also  yield  a  relative  amplitude  change  which  might  be 
beneficial.  However,  with  the  possible  exception  of  polarization  diver¬ 
sity,  the  amplitude  sensitivity  to  the  parameter  changes  is  probably 
too  small  to  be  useful.  Diversity  techniques  will  be  discussed  in 
more  detail  in  the  following  chapters. 


III. 


SUMMARY  OF  PRIOR  GLINT  INVESTIGATIONS 


A.  Early  History 

Target  scintillation  effects  were  first  observed  in  the  form  of 
target  fading  during  early  search  radar  experiments.  This  scintilla¬ 
tion,  which  was  basically  amplitude  noise,  effectively  reduced  the 
maximum  detection  range  of  the  search  radar  equipment. 

Scintillation  was  soon  observed  in  the  testing  of  experimental 
tracking  radars.  Throughout  most  of  the  World  War  II  years,  the 
effects  of  target  noise  on  tracking  performance  were  thought  to  be 
easily  explained  as  a  direct  result  of  amplitude  scintillation  intro¬ 
ducing  an  additional  modualtion  component  at  the  scanning  or  lobing 
frequency.  Such  modulation  interfered  with  the  desired  tracking 
modulation  produced  by  the  conical-scan  or  lobe-switching  mechanism. 

This  theory  was  supported  by  the  dramatic  effects  of  propeller  modu¬ 
lation  on  radar  tracking  noise  [2].  Since  this  additional  modulation 
due  to  amplitude  scintillation  affected  the  angle  tracking  loop  via 
the  scanning  frequency,  and  because  amplitude  scintillation  demonstrates 
a  decreasing  power  density  with  increasing  frequency,  it  was  logical 
to  attempt  to  reduce  the  effect  by  increasing  the  lobing  or  scanning 
rate.  In  fact,  it  was  assumed  that  if  a  tracking  mechanism  could  be 
devised  which  would  derive  tracking  information  from  all  lobes  simul¬ 
taneously  (i._e. ,  an  infinite  scanning  rate),  ic  would  be  possible  to 
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from  various  parts  of  a  complex  target.  These  time-varying  phase 
angles,  in  the  phasor  summation  performed  by  the  antenna  aperture, 
resulted  in  a  time  variation  in  the  apparent  angle-of-arrival  of  the 
radar  return.  As  early  as  1951,  Brockner  [4]  published  a  paper  dis¬ 
cussing  the  relative  relationship  of  the  various  forms  of  tracking 
noise,  jL.jJ.  amplitude  noise,  glint  noise,  servo  noise,  and  receiver 
noise.  He  concluded  that  at  relatively  short  ranges,  the  glint  noise 
was  the  predominant  term. 

While  there  was  basic  agreement  as  to  the  origin  of  glint,  there 
evolved  two  schools  of  thought  as  to  the  appropriate  modelling  technique 
and  interpretation  of  results.  These  might  be  loosely  classified  as 
the  statistical  approach  and  the  deterministic  technique. 

Statistical  methods 

Statistical  techniques  were  first  developed  by  De.’ano  (at  Hughes, 
1950)  and  later  published  in  a  paper  which  has  served  as  the  basis  for 
many  subsequent  glint  analyses  [5].  Using  as  a  model  a  set  of  iso¬ 
tropic  secondary  radiation  sources  located  near  the  tracking  axis  and 
having  statistically  independent  amplitudes  and  phases,  he  obtained 
an  expression  for  the  probability  density  of  the  fluctuation  of  the 
output  of  an  open-loop  tracker.  He  assumed  that  gross  changes  in  aspect 
wera  sufficiently  slow  to  guarantee  stationarity  of  the  assumed  statis¬ 
tics  of  the  model,  and  that  the  number  of  random  scatterers  was  large 
enough  to  produce  a  Rayleigh  distribution  of  the  echo  amplitudes. 

(Kerr  [6]  has  determined  that  as  few  as  five  or  six  sources  having 
independent  phases  and  constant  amplitudes  is  sufficient  for  Rayleigh 


Muchmore  presented  a  review  of  the  available  experimental  data  on  both 


Deterministic  Methods 

Concurrently  with  the  work  using  the  statistical  approach  to 
target  modelling,  other  investigators  were  developing  deterministic 
techniques.  At  the  University  of  Texas,  Feagin  and  Watson  [9,10]  used 
a  multipoint  target  to  obtain  the  balance  point  (the  point  the  radar 
would  track  if  its  servos  were  infinitely  fast)  for  both  a  conical- 
scan  radar  and  for  a  monopulse  radar  which  used  amplitude  in  one  track 
ing  plane  and  phase  in  the  other.  At  Ohio  State,  Peters  and  Weimer 
[11,12]  used  the  coherent  summation  of  contributions  from  a  similar 


model  and  derived  general  expressions  for  the  balance  point  of  various 


types  of  radars.  Defining  the  pointing  error  as  the  difference  between 
the  actual  balance  position  of  an  ideal  tracker  and  the  balance  position 
which  would  exist  if  the  signals  from  all  the  point  scatterers  arrived 
in  phase,  they  established  bounds  for  this  quantity  in  terms  of  the 
ratio  of  the  total  amplitudes  received  under  these  two  conditions.  The 
Ohio  State  group  also  made  extensive  measurements  of  balance  point  and 
cross-section  for  various  geometrical  shapes. 

Lewis,  [13],  further  developed  the  technique  of  Peters  and  Weimer 
to  show  that  conical-scan  and  phase  monopulse  radars  (ideally)  track 
the  phase  center  of  the  target.  This  result  was  determined  from  a 
mathematical  development  for  the  balance  point,  as  measured  from  some 
reference  point  on  the  target,  in  terms  of  the  total  field  phase  at 
the  radar.  The  balance  point  was  seen  to  vary  directly  as  the  deriva¬ 
tive  of  the  phase  with  respect  to  change  of  aspect  angle. 

While  Lewis  did  not  state  it  as  such,  the  derivative  of  phase  with 
aspect  angle  is  actually  the  slope  of  the  phase  front,  or  line  of  con¬ 
stant  phase.  The  concept  of  visualizing  angular  scintillation  as  dis¬ 
tortion  of  the  echo  phase  front  was  reported  by  Howard  of  NRL  in  late 
1959.  This  concept  has  probably  contributed  more  to  the  understanding 
of  glint  than  any  other  work.  In  his  report  [14],  Howard  showed  that 
the  slope  of  the  echo  signal  phase  front  is  identical  to  the  angular 
error  for  both  dual-  and  multi-point  targets.  He  further  reasoned  that 
since  the  distortion  is  a  property  of  the  backscattered  wave,  glint 
affects  all  angular-locating  devices,  even  search  radar. 


Analog  Simulation 


Because  it  was  evident  that  much  of  the  discrepancy  between  various 
experimental  measurements  was  due  to  the  difficulty  of  controlling  the  ex¬ 
perimental  variables,  particularly  of  an  aircraft  in  flight,  consider¬ 
able  effort  was  devoted  to  the  development  of  laboratory  glint  simula¬ 
tors.  MIT  devised  a  simulator  for  the  case  of  a  two-point-source  target 
and  an  interferometer-type  tracker  [15].  By  1956,  NRL  had  constructed 
a  closed-loop  simulator  capable  of  generating  tracking  data  for  various 
multiple-reflector  rigid  body  targets  having  a  variety  of  random  tar¬ 
get  motions  [2] .  This  simulator  was  particularly  useful  in  subsequent 
studies  of  the  effects  of  servo  bandwidth,  target  noise  bandwidth,  and 
AGC  characteristics  on  closed-loop  tracking  performance.  However,  it 
was  partially  limited  in  that  it  had  no  provision  for  non-random  rela¬ 
tive  motion  of  different  parts  of  the  target,  as  might  be  caused,  for 
example,  by  flexing  of  control  surfaces. 

AGC  Studies 

Because  the  effect  of  target  glint  is  to  introduce  a  component 
of  target  noise  into  the  tracking  loop,  several  theoretical,  simulation, 
and  experimental  investigations  were  undertaken  to  determine  the 
effect  of  automatic  gain  control  time-constant  on  this  noise  [2,7,16, 

17].  It  was  generally  agreed  that  a  slow  AGC  reduced  the  angle 
fluctuations.  However,  Howard  [17],  showed  that  a  slow  AGC  enhances 
a  new  noise  component,  proportional  to  true  tracking  errors,  which  is 
due  to  amplitude  noise  modulation  of  tracking-error  signals.  He  con¬ 
cluded  that,  because  of  the  relative  magnitudes  of  these  two  noise 


components  under  most  tracking  conditions,  a  fast  AGC  would  minimize 
over-all  tracking  noise. 

B.  Recent  Contributions 

The  first  literature  of  this  decade  appeared  in  the  form  of  a 
debate  between  the  major  advocates  of  the  two  differing  approaches 
to  glint  modelling.  Commenting  on  the  1960  unclassified  version  of 
Muchmore's  scintillation  spectra  work  [7],  Peters  and  Weimer  question¬ 
ed  the  validity  of  the  assumption  of  using  a  random  distribution  of 
scatterers  as  an  aircraft  model  [18].  They  stated  that  if  the  range 
of  aspects  is  restricted,  or  if  the  time  of  observation  is  held  to  the 
order  of  one  second  (a  typical  servo  response  time).,  the  amplitude  of 
the  echo  is  not  likely  to  be  Rayleigh  distributed.  Because  Delano 
reasoned  that  the  Rayleigh  statistics  of  observed  data  (which  was  ob¬ 
tained  over  wide  aspects)  indicated  that  five  or  six  scatterers  with 
equal  amplitudes  and  random  phases  was  an  appropriate  model,  the 
absence  of  Rayleigh  statistics  (in  an  actual  tracking  situation)  might 
invalidate  such  a  model.  Peters  and  Weimer  further  noted  that  their 
studies  of  cross-section  patterns  of  two  different  aircraft  indicated 
that  the  depth  of  successive  minima  and  the  magnitude  of  successive 
maxima  tend  to  repeat.  Such  a  repetition  would  not  exist  if  there  were 
more  than  two  effective  scatterers.  They  concluded  their  criticism  of 
the  indiscriminate  application  of  Delano's  model  by  attempting  to  demon¬ 
strate  discrepancies  which  arise  when  the  model  is  used  in  the  analysis 
of  the  two-point  target  and  linear  scatterer  cases.  In  essence,  they 
showed  that,  using  Delano's  results,  there  was  no  pointing  error  for 


the  case  of  equal  scattering  amplitudes,  except  at  null  points,  where 
the  amplitude  of  the  return  is  zero.  Because  error  was  known  to  exist, 
it  was  argued  that  this  descrepancy  illustrated  the  inadequacy  of  the 
statistical  model. 

In  the  same  issue  of  the  IRE  Transactions,  Muchmore  replied  to  the 
comments  of  Peters  and  Weimer  [19].  He  refuted  their  descrepancy  argu¬ 
ment  by  pointing  out  that,  for  the  two  source  target,  the  region  of  the 
nulls  is  exactly  the  region  where  the  glint  error  occurs.  (This  is  in 
agreement  with  the  findings  of  Howard  [14]).  Because  in  practice  the 
probability  of  all  scatterer  amplitudes  being  exactly  the  same  is 
essentially  zero,  he  stated  that  the  case  of  similar,  but  not  identical, 
amplitudes  is  more  instructive.  In  this  case,  the  nulls  do  not  go  com¬ 
pletely  to  zero,  thus  the  return  can  be  detected  and  glint  observed. 

With  regard  to  the  question  of  the  number  of  scattering  centers  effec¬ 
tive  over  small  aspect  angles,  Muchmore Ts  reply  was  more  qualitative. 

He  simply  pointed  out  that  the  assumption  of  many  random  scatterers 
was  the  most  reasonable. 

The  debate  was  continued  in  a  subsequent  issue  of  the  Transactions, 
with  communications  from  Delano,  Peters  and  Weimer,  and  Siegel  (of  the 
University  of  Michigan)  [20,21,22,23].  No  new  information  was  revealed 
by  these  discussions,  with  one  possible  exception.  Delano  noted  that 
aircraft  model  measurements  by  Radiation,  Inc.  indicated  that  echo  ampli¬ 
tudes  were  Rayleigh  distributed  over  as  little  as  5°  of  aspect  angle. 

The  importance  of  this  fact  was,  of  course,  challenged  by  Peters  and 
Weimer,  who  pointed  out  that,  because  of  typical  aircraft  velocities 


and  tracking  servo  response  times,  0.5°  was  a  more  appropriate  observa¬ 
tion  interval.  Their  contention  that  the  Rayleigh  distribution  may  not 
hold  over  such  a  small  angular  increment  is  supported  by  Katzin  [24], 
whose  presentation  of  measurements  on  a  B45  aircraft  show  that  Rayleigh- 
distributed  envelopes  are  observed  only  if  the  aspect  angle  changes  by 
more  than  four  degrees  over  the  period  analyzed.  In  an  Australian  study. 
Holt  [25]  also  observed  a  departure  from  the  Rayleigh  distribution 
during  portions  of  experimental  runs. 

There  have  been  several  glint-oriented  papers  published  in  the  open 
Russian  literature.  Ostrovityanov  [26]  observed  that  the  results  of 
several  American  investigators  had  indicated  a  possibility  that  the 
angular  error  could  become  infinite.  Because  this  is  not  a  practical 
result,  he  outlined  a  mathematical  development  which  seems  to  indicate 
that  Meade's  result  for  the  boresight  angle  [1]  should  be  instead  for 
the  tangent  of  that  angle.  (This  is  a  plausible  result,  for  the  infinite 
argument  which  arises  in  the  case  of  two  equal-amplitude,  oppositely- 
phased  sources  would  then  lead  to  a  boresight  angle  of  90° ,  in  agreement 
with  the  phase-front  concept  of  Howard  [14]).  He  also  presented  a 
detailed  statistical  development  in  which  he  refined  the  proDability 
density  expressions  of  Delano  to  account  for  certain  physical  discrep¬ 
ancies.  Gubonin  [27]  extended  Delano's  result:  to  include  the  probabi¬ 
lity  density  function  for  glint  when  there  is  one  sinusoidal,  determin¬ 
istic  reflector  included  among  the  random,  independent  reflectors. 

For  this  case,  the  density  function  is  more  narrow  than  that  for  random 
reflectors  alone,  i.e.,  it  is  less  probable  that  the  apparent  target 
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position  lies  outside  the  extent  of  the  target.  It  was  the  existence 
of  such  a  dominant  deterministic  reflector  which  Moll  believed  respon¬ 
sible  for  the  observed  departure  from  Rayleigh  statistics  [25]. 

In  1967,  the  Battelle  Memorial  Institute  began  the  task  of  perform¬ 
ing  a  complete  study  and  analysis  of  the  glint  literature  (both  clas¬ 
sified  and  unclassified)  as  applied  to  air  defence  missile  system 
design.  Their  report,  an  excellent  analysis  of  all  glint  investigations, 
includes  an  unclassified  description  of  a  "universal  glint  model"  [28], 
which  was  apparently  developed  as  a  continuation  of  the  work  of  Peters, 
et.  al. ,  at  Ohio  State. 

The  Battelle  approach  was  a  deterministic  one,  but  it  was  unique 
in  that  it  did  not  require  isotropic  sources.  Instead,  it  was  based 
upon  the  coherent  summation  of  the  fields  due  to  "scattering  complexes," 
or  extended-target  reflectors,  which  may  represent,  for  example,  air¬ 
craft  control  surfaces  or  jet  engine  intakes.  Assuming  such  complexes 
to  be  non-interacting,  the  non-isotopic  amplitude  and  phase  patterns 
of  eacn  complex  are  measured  on  a  coherent,  static  cross-section  pattern 
range,  and  the  total  target  field  obtained  by  the  summation  of  the 
appropriate  contributions  from  the  complexes.  For  cases  where  complexes 
are  interacting,  due  to  coupling  or  multiple  reflections,  these  com¬ 
plexes  are  lumped  together  and  considered  as  a  single  complex.  Using 
the  summation  expression  for  the  total  field,  Battelle  obtained  the 
direction  of  the  normal  to  the  phase  front  by  taking  the  gradient  of 
the  phase  term.  The  spatial  derivatives  of  the  amplitude  and  phase  com¬ 
ponents  of  the  individual  complexes,  which  result  from  the  gradient 
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operation,  are  obtained  from  the  measured  cross-section  plots.  Relative 
motion  between  complexes  is  included  by  adjusting  the  individual  phase 
terms  to  include  differential  Doppler. 

In  discussing  the  model,  Battelle  reached  one  conclusion  which 
differed  from  a  generally-accepted  hypothesis  concerning  the  correla¬ 
tion  between  glint  errors  and  echo-amplitude  minima.  Howard  [17]  and 
others  had  postulated  that  there  was  a  negative  correlation  between 
these  quantities,  jL.e^  glint  error  was  maximum  at  regions  of  local 
minima.  This  correlation  was  based  upon  an  extension  of  the  two-point- 
source  case.  The  Battelle  model  indicated  that  a  field  minima  from 
any  complex  would  result  in  large  glint  error.  Such  a  component  minima 
would  not,  in  general,  correspond  to  a  general  amplitude  null. 

The  most  recert  contribution  to  the  glint  literature  was  from 
Dunn  and  Howard  of  NRL  [29].  Their  objective  was  to  show  that  glint 
for  a  complex  target  corresponds  to  a  deviation  of  the  Poynting  vector 
from  the  radial  direction,  and  thus  a  corresponding  deviation  in  the 
direction  of  the  power  flow.  To  demonstrate  this,  they  assumed  a 
target  model  consisting  of  a  collinear,  non-uniform  array  of  electrically 
short  dipoles.  Because  the  spatial  patterns  of  the  electric  and  magne¬ 
tic  fields  from  such  an  element  are  well  known,  it  was  mathematically 
feasible  to  determine  the  total  E  and  H  fields  resulting  from  the  coherent 
summation  of  the  individual  element  fields  at  a  general  observation 
point.  The  analysis  using  these  total  fields  showed  that  there  was  indeed 
a  component  of  the  time-averaged  Poynting  vector  orthogonal  to  the 
radial  direction.  This  proved  that  the  power  flow  is  not  necessarily 
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radial,  and,  at  a  given  observation  point,  may  be  in  such  a  direction 
as  to  indicate  the  target  position  to  be  outside  the  actual  target  span 
Because  of  this,  even  the  simple  search  radar  is  susceptible  to  glint. 
This  was  illustrated  by  observing  the  change  in  the  measured  pattern 
of  a  horn  antenna,  as  the  relative  phase  between  a  closely-spaced  pair 
of  illuminating  horns  was  varied.  The  authors  discussed  the  glint 
averaging  effect  of  a  finite  antenna  aperture,  and  developed  equations 
for  relationships  between  Doppler  and  angular  scintillation.  They  also 
discussed  the  effects  of  various  diversity  techniques;  these  comments 
will  be  further  reported  in  the  next  section. 

C.  Application  of  Diversity 

As  early  as  1956,  Peters  and  Weimer  [11]  suggested  two  diversity 
techniques  which  might  be  used  to  reduce  pointing  errors.  Both  of 
these  seem  to  be  based  on  the  concept  that  maximum  error  occurs  at 
amplitude  minima.  The  first  suggestion  was  to  servo-control  the  radar 
frequency  to  maintain  the  amplitude  of  the  return  within  certain  limits 
The  other  idea  was  to  accept  and  clamp  the  error  signal  to  the  servo 
only  when  the  amplitude  is  greater  than  a  certain  minimum  value. 

Birkemeier  and  Wallace  [30]  investigated  the  case  of  pulse-to- 
pulse  frequency  modulation  for  possible  angle  tracking  accuracy  improve 
ment.  They  used  a  statistical  target  model  consisting  of  elements 
uniformly  distributed  along  a  strip  and  having  independent  random  phase 
angles.  Their  motivation  in  performing  the  study  was  that  the  glint- 
inducing  phase  interference  effects  from  the  target  could  be  made  to 
occur  more  rapidly  by  progressively  varying  the  radar  carrier  frequency 


Ray  [32]  analyzed  the  improvement  in  radar  range  and  angle  detection 
with  frequency  agility.  He  stated  that  the  amount  of  frequency  shift 
necessary  for  decorrelation  of  the  returns  from  a  target  is  a  function 
of  the  range  dimensions  of  the  target  complex.  He  analyzed  both  the  case 
where  the  target  range  dimension  is  much  less  than  the  range  resolution 
of  the  radar,  and  the  case  where  this  dimension  approaches  or  exceeds 
the  range  resolution.  In  both  cases,  a  significant  improvement  in 
performance  was  predicted.  While  his  discussion  of  angle  detection  was 
restricted  to  the  angular-locating  accuracy  of  a  search-type  beam,  much 
of  his  analysis  should  be  applicable  to  automatic  tracking  radars. 

The  most  pertinent  discussion  of  diversity  (as  applied  to  improved 
radar  tracking)  to  appear  in  the  open  literature  was  contained  in  the 
recent  paper  of  Dunn  and  Howard  [29].  The  authors  first  commented  on 
space  diversity,  as  achieved  by  locating  several  radar  receivers  at 
positions  widely  separated  with  respect  to  the  phase-front  distortion 
pattern.  Each  receiver  then  observes  independent  samples  of  angular  scin¬ 
tillation.  The  average  of  the  various  error-signal  outputs  thus  ex¬ 
hibits  a  reduction  of  target  noise  power.  The  authors  point  out  that 
the  averaging  provided  by  an  antenna  when  the  phase-front  distortion  region 
spans  a  large  portion  of  the  antenna  aperture  is  also  a  form  of  space 
diversity.  Frequency  diversity  is  indicated  as  another  technique  for 
averaging  target  scintillation,  either  instantaneously  or  within  a  short 
period  of  time.  The  object  is  to  provide  sufficient  separation  between 
frequencies  such  that  the  scintillation  at  each  frequency  is  uncorrelated 
with  that  at  the  other  frequencies.  An  example  presented  shows  that, 
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for  two  reflectors  displaced  by  10  yards  in  range,  a  frequency  change 
of  18  MHz  causes  a  full  cycle  of  change  in  the  received  relative 
phase.  The  authors  note  that  space-diversity  effectiveness  is  determined 
by  the  angle  subtended  at  the  target,  while  frequency-diversity  effective¬ 
ness  is  determined  by  the  range  distribution  of  the  target.  It  may 
therefore  be  possible  to  combine  effectively  the  two  techniques.  In  the 
case  of  polarization  diversity,  the  effectiveness  is  dependent  upon  the 
polarization  characteristics  of  each  scattering  element.  It  was  also 
noted  that  some  of  the  deeper  signal  fades  may  be  averaged  by  diversity 
techniques . 


D.  Elevation  Error  Studies 

There  exists  another  problem  in  radar  tracking  which  has  an  effect 
quite  similar  to  that  of  target  glint.  This  problem  is  the  multipath 
phenomenon,  which  produces  elevation  errors  when  a  radar  is  tracking  a 
target  at  a  low  elevation  angle  above  the  earth.  Because  the  radar 
receives  a  return  from  the  target  image  in  the  (partially)  conducting 
earth's  surface,  in  addition  to  the  direct  return  from  the  target,  and 
because  these  returns  are  from  slightly  different  angles,  the  effect  is 
quite  similar  to  the  glint  caused  by  the  two-reflector  target.  There 
have  been  many  analytical  and  theoretical  studies  of  this  glint-related 
problem.  Much  work  has  been  done  in  England  by  the  RRE  and  others  [33, 

34]  and  the  effect  is  discussed  in  several  texts  and  reports  (e.£.  [6,35]). 


E.  Comments 

The  preceding  summary  of  investigations  is  not  intended  to  be  an 


exhaustive  survey  of  all  glint  studies.  However,  it  is  felt  the  material 
covers  the  major  contributions,  as  reported  in  the  open  literature,  in 
sufficient  detail  to  indicate  the  direction  of  the  glint  research  efforts. 

Several  conclusions  might  be  reached  concerning  prior  glint  studies. 
It  appears  that  the  various  investigators  have  lacked  unity  of  purpose 
in  the  design  and  techniques  of  their  analytical  and  experimental  efforts. 
This  is  no  doubt  due  to  the  extended  period  of  time  during  which  the 
studies  were  conducted,  and  the  varying  research  goals  of  the  various 
facilities  doing  the  research.  However,  in  some  cases  it  may  be  that 
the  research  was  performed  to  substantiate  the  need  for  a  particular 
radar  design  or  tracking  technique;  in  such  cases  the  common  research 
goal  of  increased  understanding  of  a  problem  may  have  suffered. 

The  lack  of  unity  of  purpose  of  the  various  researchers  is  most 
evident  in  the  reports  of  experimental  series.  The  parameter  studies 
have  ranged  from  true  glint  (property  of  the  propagating  wave  front) 
to  tracking  noise,  as  measured  by  both  open-loop  and  closed-loop  receiver 
output  voltages  and  with  and  without  amplitude  scintillation  effects 
removed.  In  addition,  the  data  may  be  referenced  to  a  fixed  point  on 
the  target  or  to  the  radar  "center  of  gravity,"  and  may  or  may  not 
include  the  effects  of  boresight  errors.  In  many  cases,  it  is  not  pos¬ 
sible  to  ascertain  which  of  these  forms  of  the  glint  parameter  was 
measured.  Such  vast  differences  make  it  extremely  difficult  to  compare 
the  various  results.  For  this  reason,  considerable  planning  should  go 
into  the  development  of  any  future  measurement  programs  to  insure  that 
the  effect  being  measured  is  glint  alone,  and  that  the  results  depend 


as  little  as  possible  upon  the  radar  used  in  the  measurements.  In  fact, 
it  might  be  more  efficient  to  postpone  extensive  test  programs  until 
more  analytical  model  computations  have  been  performed.  An  experimental 
effort  could  then  be  devoted  to  verification  of  these  computations  by 
measurements  at  a  few  carefully  selected  sample  points. 

With  regard  to  the  debate  of  statistical  V£  deterministic  modelling 
techniques,  it  is  probable  that  a  better  model  is  one  which  incorporates 
features  from  both  types.  Peters  and  Weimer  argued  that  the  requirement 
for  six  or  more  scatterers  to  be  active  over  all  aspect  angles  is  too 
stringent.  This  seems  to  be  a  valid  point.  Additionally,  the  statistical 
model  demands  stationary  statistics.  However,  Delano's  contention 
concerning  the  impracticability  of  knowing  the  exact  amplitude  and  phase 
contributions  from  all  scattering  points  on  a  moving  target  also  seems 
valid.  Howard  has  noted  that  random  roll,  pitch,  and  yaw  of  an  aircraft 
target  in  flight  can  only  be  described  statistically.  Thus,  although 
it  seems  that  Delano's  statistical  model  may  not  be  valid  for  real 
targets  (because  of  the  number-of-scatterers  and  stationarity  limitations), 
it  would  be  most  desirable  to  be  able  to  include  regular  and  random 
relative  or  rigid-body  motions  of  the  target  into  a  deterministic  model 
that  does  not  require  the  assumption  of  isotropic  point  scatterers. 

The  technique  of  coherent  summation  of  the  contributions  from  scattering 
complexes,  as  introduced  by-Battelle,  provides  a  basis  for  such  a  model. 
Because  the  contributions  of  the  various  complexes  are  measured  quantities, 
the  computed  return-field  configuration  is  realistic.  The  statistical 
nature  of  the  field  when  the  target  has  relative  or  rigid-body  motion 


might  then  be  determined  analytically  be  appropriately  modifying  the 
phase  terms  of  the  individual  contributions.  Xu  this  way,  no  statistical 
stationarity  is  required,  and  the  basic  non-isotropic  nature  of  the 
scatterers  is  preserved. 

The  remainder  of  this  study  will  be  devoted  to  the  development 
and  implementation  of  an  analytical  glint  model  which  will  be  based  upon 
the  Battelle  concept  of  coherent  summation  of  the  contributions  from 
non- isotropic  scattering  complexes.  The  method  of  obtaining  the  glint 
information  from  the  resultant  total  field,  however,  will  be  considerably 
more  straightforward  than  the  Battelle  gradient  approach,  and  one  which 
is  easily  implemented  on  a  digital  computer. 


IV.  DEVELOPMENT  OE  A  GLINT  MODEL 


It  was  noted  near  the  end  of  Chapter  III  that  the  Battelle  concept 
of  a  target  model  composed  of  independent,  nonisotropic  scattering 
complexes  seems  to  be  a  very  realistic  one.  However,  the  mathematical 
processing  of  the  measured  scattering  complex  data,  as  developed  by 
Battelle,  involves  the  gradient  of  the  phase  term.  It  thus  requires 
not  only  the  measured  amplitude  and  phase  patterns  of  the  complexes, 
but  the  slope  of  these  patterns  in  orthogonal  directions  at  each  point 
in  space  as  well.  The  resulting  expression  for  the  direction  of  the 
normal  to  the  phase  front  is  rather  complex.  Moreover,  the  actual 
amplitude  and  phase  of  the  total  field  is  not  a  direct  result  of  the 
process  of  evaluating  this  expression.  It  is  reasonable  to  presume 
that  a  deterministic  model  which  would  characterize  the  amplitude 
and  phase  of  the  reflected  wave,  in  addition  to  the  normal  to  the 
phase  front,  might  prove  very  useful  in  formulating  a  signal  processing 
scheme  to  reduce  glint  error.  This  can  be  achieved  by  forming  the 
coherent  sum  of  the  contributions  from  the  scattering  complexes  (with 
the  phases  modified  for  motion,  if  required)  and  then  using  simple 
geometrical  relationships  to  determine  the  glint  error.  Information 
concerning  dynamic  effects,  such  as  the  spectra  of  amplitude  and 
angle  scintillations  during  a  particular  type  of  target  or  tracker 
maneuver,  can  be  determined  by  an  appropriate  point-by-point  analysis. 


It  should  be  evident  that  an  exact  representation  of  the  glint 
signature  at  various  aspect  angles  for  a  single  target  will  have  little 
resemblance  to  the  exact  values  at  the  same  angles  for  another  target, 
even  if  the  two  are  quite  similar.  The  precise  structure  of  the  back- 
scattered  interference  pattern  is  too  dependent  upon  target  surface 
configuration  to  permit  such  a  simple  extrapolation  of  results.  Indeed 
because  of  the  extremely  fine  lobe  structure  of  the  cross-section 
pattern  for  a  large,  complex  target,  the  amount  of  data  and  the  number 
of  calculations  necessary  to  form  a  complete  sphere  of  information 
even  for  a  single  such  target  may  be  prohibitive.  Also,  the  radar 
wavelength  at  microwave  frequencies  is  very  small  compared  to  typical 
distances  between  scattering  complexes.  An  actual  target  in  flight 
probably  has  flexing  motion  between  complexes  amounting  to  appreciable 
portions  of  a  wavelength.  Thus  the  task  of  determining  the  exact 
coordinates  of  the  complexes  with  respect  to  the  overall  airframe  of 
a  particular  target  cannot  be  accomplished. 

It  might  seem  that  the  limitations  just  discussed  would  negate 
any  utility  of  a  deterministic  model.  However,  it  is  emphasized  that 
the  motivation  in  the  development  of  the  model  to  be  described  is  to 
permit  the  exact  computation  of  representative  field  parameters  over 
representative  aspect  sectors  for  representative  classes  of  targets. 
Then  a  determination  of  the  dynamic  behavior  of  these  parameters 
along  a  realistic  scries  of  sample  points  within  each  sector  can  be 
made.  Such  a  procedure  will  permit  the  analytical  evaluation  of 
various  compensation  and  smoothing  schemes  in  a  set  of  realistic  field 


environments . 


Several  basic  assumptions  are  necessary  before  beginning  the 
development  of  the  model.  It  is  presumed  that  the  technique  of 
representing  a  radar  target  by  an  array  of  scattering  complexes  is 
valid.  This  assumption  demands  that  the  scattering  pattern  of  each 
complex  be  independent  of  other  complexes,  and  that  insignificant 
error  be  introduced  by  the  physical  dissection  of  the  target  into 
scattering  complexes,  £.&.  edge  effects  at  lines  of  dissection. 
Moreover,  it  is  necessary  that,  at  ranges  of  interest  for  the  complete 
target,  the  fields  of  the  individual  complexes  have  a  far-field 
amplitude  and  phase  variation  with  range.  (The  validity  of  this 
assumption  will  be  investigated  in  a  following  section.)  The 
discussion  will  be  oriented  toward  the  problem  of  a  monostatic 
tracker  in  a  homing  situation;  however,  the  techniques  can  be  applied 
to  any  tracking  radar,  although  a  bistatic  configuration  would  be 
more  involved. 


A.  Obtaining  Amplitude  and  Phase  Patterns 
for  the  Scattering  Complexes 

For  maximum  versatility  of  the  model  in  studying  the  effects  of 
diversity  on  glint  behavior,  it  is  desirable  that  the  amplitude  and 
phase  data  for  individual  scattering  complexes  be  recorded  with 
frequency  and  polarization  as  variable  parameters.  (A  discussion  of 
space  diversity  will  be  delayed  until  a  later  section.)  A  possible 
scattering  range  implementation  for  obtaining  the  data  is  depicted 
schematically  in  Fig.  4-1.  This  configuration  permits  the  recording 


of  the  amplitude  and  phase  of  orthogonal  components  of  the  reflected 
signal,  as  a  function  of  the  aspect  angle  of  the  scattering  complex 
and  the  transmitted  polarization  and  frequency. 

The  transmitter  comprises  a  backward  wave  oscillator  (BWO)  which 
is  frequency -stepped  by  an  analog  voltage  from  a  frequency  programmer. 
For  each  frequency  increment,  a  polarization  programmer  alternately 
connects  the  BWO  to  horizontally  and  vertically  polarized  feeds  on  the 
transmitting  antenna  by  means  of  a  SFDT  diode  switch  (SW) . 

Horizontally  and  vertically  polarized  feeds  on  the  receiving 
antenna  provide  simultaneous  signals  to  separate  phase-and-amplitude 
receivers.  (A  single  receiver  could  be  used  on  a  time-shared  basis.) 

A  sample  of  the  BWO  output  provides  the  required  phase  reference  for 
the  receivers. 

The  target  positioner  has  two  degrees  of  freedom  about  some 
defined  coordinate  origin  on  the  scattering  complex  (azimuth  and 
elevation  angle  variations)  and  is  controlled  by  a  positioner  pro¬ 
grammer.  This  programmer,  together  with  the  frequency  and  polarization 
programmers,  provide  digitally-coded  outputs  identifying  the  aspect 
angle,  frequency,  and  polarization  parameters.  This  data  is  recorded 
on  magnetic  tape,  along  with  the  digitized  amplitude  and  phase  outputs 
of  the  horizontal  and  vertical  receivers .  Synchronism  of  the  various 
parameter  changes  is  maintained  by  a  clock. 

Several  comments  regarding  the  range  implementation  just  discussed 
are  in  order.  It  is  envisioned  that  the  coherent  static  measurements 
would  be  made  on  a  large,  ground-level  radar  range  such  as  the  RAT  SCAT 
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facility  in  New  Mexico  [36],  Because  such  a  range  utilizes  the 
multipath  effect . in  achieving  target  illumination  and  receiving  the 
radar  return,  it  actually  has  a  finite  bandwidth  associated  with  a 
particular  combination  of  antenna  and  target  heights.  Since  the 
proposed  range  tests  involve  measurements  over  a  band  of  frequencies, 
a  careful  study  should  be  made  to  verify  that  this  band  is  within 
the  operating  bandwidth  of  the  range. 

Because  the  validity  of  this  model  depends  upon  the  summation 
of  appropriate  field  contributions  from  the  various  complexes,  the 
registration  accuracy  with  which  each  complex  is  mounted  on  the 
positioner  pylon  will  be  of  extreme  importance.  Also,  care  should 
be  taken  to  be  sure  that  the  physical  size  of  the  complex  is  within 
the  limits  imposed  by  the  range  length  and  antenna  diameter.  If  the 
data  is  to  be  essentially  independent  of  the  type  of  antenna  used  in 
the  measurements,  uniform  illumination  across  the  aperture  is  desired. 
Caster  137]  states  that  if  a  complex  target  subtends  less  than  1/8 
beamwidth,  this  condition  will  be  achieved.  Thus,  using  the  approximate 
relationship  for  beamwidth  of  a  uniformly  illuminated  aperture, 

0^  =  A/D,  a  1  meter  dish  operating  at  10  GHz  would  permit  a  maximum 
target  dimension  of  3.75  meters  at  a  range  of  1  km. 

Transmitting  alternate  orthogonal  linear  polarizations,  as  opposed 
to  continuous  circular  polarization,  is  preferred  because,  with  this 
measurement  technique,  the  phase,  of  the  received  signal  is  not  lost, 
as  in  the  case  of  a  conventional  range  receiver  using  a  bolometer 
(integrating)  detector.  The  rotation  frequency  of  the  linear  electric 
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rotation  of  the  two  degrees  of  freedom  for  the  scattering  range  posi¬ 
tioner.  Assume  the  range  to  the  positioner  to  be  constant  for  all 
scattering  measurements  and  let  frequency  be  a  fixed  parameter.  Then, 
for  a  specified  receiving  polarization,  the  amplitude  and  phase  data 
are  a  function  of  the  coordinates  6*  and  <£■*■,  as  measured  in  the  system 
centered  at  0^.  Now  let  a  spherical  coordinate  system,  with  origin  0, 
be  associated  with  some  point  on  the  target.  The  aspect  point  F,  at 
coordinates  (r,0,4>),  and  thus  all  field  values  determined  from  the  model, 
will  be  referenced  to  this  origin.  Finally,  let  the  location  of  each 
origin  0i,  relative  to  0,  be  specified  by  the  coordinates  (d^.e^,^) , 
and  let  all  coordinate  axes  have  the  same  relative  orientation.  This 
geometry  is  shown  in  Fig.  4-2. 

Propagation  path  determination 

Because  the  technique  of  coherent  summation  uses  the  amplitude  and 
phase  contribution  from  each  complex  at  the  field  point  P,  it  is  neces¬ 
sary  to  know  the  spatial  attenuation  and  phase  retardation  over  the 
propagating  path  from  each  complex.  These  effects  can  be  easily  deter¬ 
mined  if  the  distance  from  P  to  each  secondary  origin  is  known.  To 
derive  an  expression  for  this  distance,  consider  the  geometry  shown  in 
Fig.  4-3,  where  the  secondary  origin  is  denoted  by  0".  Ftv-m  the  law 
of  cosines, 

d2  =  r2  +  d2  _  2rd  cos  a  (4-1) 

or 

D-^  =  [r^  +  d^  -  2rd  cos  . 

=  (1/ r)  [  1  +  (d/r)2  -  2(d/r)  cos  al^--^) 


(4-2) 
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Letting  x  =  cos  a  and  G  =  d/r, 

D-1  =  (1/r) [1  +  G2  -  2Gx] .  (4-3) 

This  can  be  expanded  in  a  binomial  series,  provided 

(G2  -  2Gx)2  <  1.  (4-4) 

To  express  the  convergence  condition  in  terms  of  physical  parameters, 
note  that,  from  (4-3), 

D2  =  r2  +  r2(G2  -  2Gx) 
or 

G2  -  2Gx  -  (D/r)2  -  1.  (4-5) 

Thus  (4-4)  is  equivalent  to 

[(D/r)2  -  l]2  <  1 
or 

(D/r)4  -  2 (D/r) 2  <  0 
or 

(D/r)  <  Ji  .  (4-6) 

Inequality  (4-6)  gives  the  convergence  condition  in  terms  of  the  distances 
from  the  primary  and  secondary  origins  to  the  field  point  P.  The  quantity 
(D/r)  is  positive,  and  for  a  given  r  and  d,  has  a  maximum  value  when  0" 
and  P  are  diametrically  opposed,  or 

Dmax  "  r  +  d.  (4-7) 


Substituting  this  maximum  value  in  (4-6),  the  maximum  displacement  of 


Fig.  4-3. — Geometry  for  deriving  the  path  length  D  and  for 
determining  the  coordinates  of  P  relative  to  the 
secondary  origin  0''. 
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the  secondary  origin  from  the  primary  origin  becomes 

W  =  0.414r.  (4-8) 

This  condition  is  easily  satisfied  in  practice  because  the  range  is 
usually  much  greater  than  the  maximum  target  dimension.  Returning  to 
(4-3) ,  the  binomial  series  expansion  is 

D-1  =  (l/r)[l-(l/2)(G2-2Gx)+(3/8)(G4-4G3x2-8G3x3) 

-  (15/48) (G6-6G5x+12G4x2-8G3x3)-  ...].  (4-9) 

Collecting  like  powers  of  G, 

D”1  =*  (l/r)  [1  +  Gx  +  G2(  1/2)  (3x2  -  1) 

+  G3(1/2)(5x3  -  3x)  +  ...].  (4-10) 

Replacing  G  by  d/r,  (4-10)  is  seen  to  be  of  the  form 

00 

D"1  =  (l/r)  l  (d/r)£  PA(x)  (4-11) 

£=0 

where  P&(x)  is  the  Legendre  polynomial  of  order  i  and  argument  x. 

In  order  to  use  the  expression  (4-11)  to  find  the  distance  D,  it 
is  necessary  to  determine  the  argument  x  =  cos  a  in  terms  of  the  spheri¬ 
cal  coordinates  of  P  and  0".  In  Fig.  4-3,  consider  the  radial  lines 
to  0"  and  P  to  be  the  vectors  d  and  r,  respectively.  Forming  the  scalar 
product  of  these  vectors, 

r  *  d  =  rd  cos  a  =  xx'  +  yy'  +  zz'  (4-12) 
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where  (x,y,z)  and  (x  ,y  ,z  )  are  the  rectangular  coordinates  of  P  and 
0",  respectively.  Substituting  the  relationships  between  spherical  and 
rectangular  coordinates, 
x  «  r  sin  0  cos  <j> 
y  =  r  sin  0  sin  <£ 

z  =  r  cos  <j>  (4-13) 

in  (4-12)  and  solving  for  cos  o  yields 

cos  a  =  cos  0  cos  0*  +  sin  0  sin  0'  cos  (<|>-<j>')  (4-14) 

where  (0 , 4>)  and  (0’,$')  are  the  spherical  coordinates  associated  with 
P  and  0",  respectively. 

Transformation  to  secondary  coordinates 

To  find  the  total  field  at  P,  it  is  necessary  to  know  the  phasor 
contribution  of  each  complex  C-^.  Each  contribution  is  found  by 
examining  the  recorded  data  which  represents  the  spatial  patternf  of 
the  complex,  extracting  from  the  data  the  amplitude  and  phase  values  in 
the  direction  given  by  the  angular  coo ruinates  of  P  relative  to  0-^,  and 
modifying  these  values  to  include  the  effects  of  spatial  attenuation 
and  phase  retardation  over  the  appropriate  propagation  path.  Because 
the  angular  coordinates  of  P  are  known  only  in  the  primary  coordinate 
system,  a  coordinate  transformation  is  required  to  obtain  the  angular 
coordinates  of  P  relative  to  Oi.* 

*The  radial  distance  to  P  is  usually  much  larger  than  the  distances 
between  the  primary  and  secondary  origins.  This  fact  may  permit  the  use 
of  the  primary  angular  coordinates  of  P  in  extracting  the  amplitude  and 
phase  values  from  the  recorded  data,  thus  avoiding  the  coordinate  trans¬ 
formations  to  each  secondary  system.  However,  an  examination  of  the  fine 
grain  structure  of  the  data  for  many  complexes  would  be  necessary  to 
validate  such  an  approximation. 


§3 


43 

Consider  again  the  geometry  shown  in  Fig.  4-3,  where  the  secondary 
origin  is  denoted  as  0".  Let  the  spherical  coordinates  of  P  with  respect 
to  0  be  (r , 8 , cp)  and  with  respect  to  0"  be  (D,8", <j>") ,  and  let  the  coordi¬ 
nates  of  0”  with  respect  to  0  be  (d>0,,<J>1).  Let  the  corresponding 
rectangular  coordinates  be  (x,y,z),  (x",y':,z")>  and  (x’jy’jZ*),  respec¬ 
tively.  The  known  parameters  are  (r,0,(j>),  (d , 0 ' , ^> * ) ,  and  D  (from  (4-11)) 
It  is  desired  to  find  the  equivalent  angular  coordinates  (0",<j>").  By 
translation  of  rectangular  coordinates, 


□ 


i 

>•* 


L 

t 

1  .‘i 


z  -  z 


X  -  X 


(4-15) 


Substituting  the  spherical- rectangular  relationships  (4-13)  on  the 
right  hand  sides  yields 

z"  =  r  cos  0  -  d  cos  8' 

x"  =  r  sin  0  cos  <j>  -  d  sin  0'  cos  <j>’.  (4-16) 


Again  using  (4-13) , 

0"  =  arccos  [z"/D] 

=  arccos  [ (r  cos  0  -  d  cos  6')/D] 

and 

<j>"  =  arccos  [x"/  (D  sin  8")  ] 


=  arccos 


r  sin  8  cos  4>  —  d  sin  81 


cos 


D  sin  0" 


t 

* 


(4-17) 


(4-18) 


Equations  (4-17)  and  (4-18)  are  the  desired  transformations. 
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Summation  of  contributions 

Let  the  amplitude  and  phase  contributions  from  complex  C^,  in  the 
direction  given  by  the  angular  coordinates  of  P  relative  to  0^,  be  Af 
and  ,  respectively.  These  are  the  values,  taken  from  the  recorded 
static  range  data,  which  represent  the  pattern  of  the  complex  in  this 
direction.  If  the  data  are  recorded  according  to  the  maximum- target- 
dimension  criterion  discussed  in  Section  A,  the  amplitude  and  phase 
will  have  a  far-field  variation  with  range.  Thus  the  amplitude  and 
phase  contributions  of  complex  at  the  field  point  P  are 

AJ  -  (4-19) 

Di 

and 

a’  =  -  2k(Di  -  RJ  (4-20) 

where  k  =  2tt / A  is  the  wave  number,  Dj.  is  the  distance  from  0^  to  P,  and 
Rfl  is  the  range  at  which  the  static  measurements  are  made  (constant 
for  all  complexes) .  The  factor  of  two  in  the  propagation  path  term 
arises  because  of  the  assumed  monostatic  configuration. 

Let  the  total  number  of  complexes  be  N>  and  let  Ei=  A| /« j.  be  the 
phasor  representation  of  the  contribution  of  Cj_  at  P.  Then  the  total 
received  signal  ut  P  is 
N  _ 

E=A/a=  l  E..  (4-21) 

i=l 

In  terms  of  known  quantities,  (4-21)  becomes 


Ih  -  2k(Dj  -•  V  . 


-  /. 

i=l 


The  procedure  for  using  (4-22)  to  find  the  total  field  E  at  a  given 
point  P  thus  involves  the  use  of  (4-11)  to  find  the  distance  to  each 
complex  C^,  the  use  of  transformations  (4-17)  and  (4-18)  to  find  the 
angular  coordinates  of  P  with  respect  co  each  origin  0-^,  and  the  use 
of  these  coordinates  to  determine  the  contributions  A^  and  from  the 
recorded  data. 

Provisions  for  motion 

Providing  for  a  turning  movement  of  the  target  requires  two  modifi¬ 
cations  of  the  above  procedure.  First,  because  the  target  is  turning, 
the  aspect  angle  is  changing  with  time.  This  change  can  be  accounted 
for  by  making  modifications  of  the  angular  coordinates  of  P  which 
reflect  the  change  in  aspect  angle.  A  derivation  of  the  coordinate 
changes  for  the  case  of  constant  roll,  pitch,  and  yaw  rates  is  given  in 
Appendix  A.  Second,  the  turning  motion  results  in  a  Doppler  shift  in 
the  radar  frequency  of  the  return  from  each  complex.  On  an  instantaneous 
basis,  such  a  frequency  change  is  equivalent  to  a  phase  change.  Thus 
the  Doppler  effect  can  be  included  an  appropriate  modification  of 
the  phase  angle  of  each  complex.  A  discussion  of  this  modification 
for  the  case  of  constant  turning  rates  is  included  in  Appendix  A. 

If  there  exists  relative  motion  between  a  complex  and  the  overall 
airframe,  as  in  the  case  of  vibratory  motion  of  a  control  surface,  the 
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phase  contribution  from  that  complex  must  be  modified.  The  appropriate 
change  consists  of  the  addition  of  a  phase  term  which  represents  the 
particular  phase  modulation  due  to  the  motion. 

C.  Determination  of  the  Glint  Error 

For  this  analysis,  the  glint  error  will  be  defined  as  the  angular 
deviation  of  the  instantaneous  apparent  radar  center  from  the  origin  0. 
Referring  to  Fig.  4-4,  the  position  of  the  apparent  radar  center,  C, 
may  be  expressed  in  linear  units  at  the  target  along  line  OC,  as 
measured  in  the  plane  JKLM,  which  is  normal  to  OP  and  contains  0.  The 
plane  OPQ  is  defined  as  the  elevation  plane,  while  the  transverse  plane 
is  that  which  is  orthogonal  to  the  elevation  plane  and  intersects  it 
along  the  line  OP.  The  glint  error  will  be  determined  separately  in 
these  two  planes  and  then  the  results  combined  to  find  the  total  error. 

The  phase  front  may  be  defined  as  the  range  change  increment,  in 
linear  units,  required  to  reach  a  position  having  some  reference  phase 
aQ.  The  phase,  a,  is  related  to  range,  r,  by 

a  =  -kr.  (4-23) 

If  the  reference  phase  is  taken  to  be  aQ  =  0,  then  the  phase  front,  as 
measured  from  the  range  coordinate  where  a  is  determined,  is  simply 

F  =  -a/k.  (4-24) 

Because  the  range  of  values  for  the  phase  angle  is 


0  <  a  <  2tt 


the  phase  front  is  a  number  which 


0  <  F  <  A.  (4-25) 

It  is  assumed  that  the  phase  front  is  a  well-behaved  function 
of  aspect  angle  in  the  vicinity  of  P.  This  will  be  said  to  be  true  if 
a  region  tihe  size  of  a  typical  airborne  radar  aperture  (less  than  one 
meter)  spans  no  more  than  one  full  cycle  of  the  distorted  wave  front. 
For  the  case  of  the  two- ref lector  target,  Dunn  and  Howard  [29] 
showed  that  this  requires  a  range  greater  than 

>w  ■  <4-26> 

where  d  is  the  reflector  spacing  and  A  is  the  aperture  diameter.  It 
is  assumed  that  this  relationship  is  also  reasonably  valid  for  more 
general  target  configurations,  and  moreover,  that  the  phase  front  is 
.  approximately  a  linear  function  of  aspect  angle  over  the  region  around 
P.  Stated  another  way,  it  is  assumed  that  the  range  is  large  enough 
for  the  phase  front  over  the  region  of  interest  to  be  approximately 
planar. 

To  determine  the  glint  error,  consider  the  geometry  depicted  in 
Fig.  4-5,  where  OC’P  lies  in  either  the  elevation  or  the  transverse 
plane,  and  where  C'  is  the  projection  of  the  apparent  radar  center  C 
in  this  plane.  The  range  R  is  assumed  to  be  much  larger  than  A  (the 
diameter  of  the  region  of  interest),  and  thus  the  contributions  and 
may  be  considered  constant  over  the  solid  angle  subtended  by  A  at 
each  origin  0^.  F-,  and  F2  are  the  values  of  the  phase  front  at  the 


two  extremes  of  A.  The  linear  glint  error  OC'  is  the  parameter  to  be 
determined. 


By  definition  of  the  apparent  target  center,  the  line  PC*  is  normal 
to  the  phase  front.  The  slope  of  the  phase  front  is 

ml  =  (F1  '  F2} M  (4"27) 

while  the  slope  of  PC’  with  respect  to  OP  is 

m2  -  OC'/R.  (4-28) 

By  similar  triangles,  these  slopes  are  identical,  therefore 

PC'  F2  -  F1 
R  A 


or 

PC'  -  R(F2  -  Fx)/A. 


(4-29) 


Substituting  from  (4-24)  for  the  phase  front  in  terms  of  the  total  field 


phases  at  the  extremes  of  the  sampling  region,  a'x  and  a!?, 

R(c i'i  -  o«) 

0C'  "  —  IT" “  ‘ 


(4-30) 


However,  in  determining  these  total  field  phases,  the  phase  contribution 
from  each  complex,  c<x,  must  be  obtained  in  a  slightly  different  manner 
than  that  indicated  by  (4-20) .  The  illumination  path  from  P  to  the 
complex  is  still  p^,  but  the  path  for  the  radar  return  from  the  complex 
is  now  the  distance  from  0X  to  the  extreme  of  the  sampling  region  where 
the  field  is  to  be  determined.  Denoting  this  distance  by  Dx, 


implementation  to  study  samples  of  this  effect  and  others  will  be 
described  in  the  next  chapter. 


V.  IMPLEMENTATION  OF  THE  MODEL 


As  an  extension  of  the  development  of  the  glint  model  described 
in  the  previous  chapter,  a  series  of  FORTRAN  IV  computer  programs 
were  written  to  illustrate  the  determination  of  target  glint  from  the 
backscattered  wavefront.  The  known  results  of  several  configurations 
of  isotropic  scatterers,  observed  in  a  far-field  situation,  were  used 
to  establish  the  validity  of  the  computer  techniques.  An  effort  was 
made  to  keep  the  programs  sufficiently  general  to  permit  their  use  in 
describing  the  effect  of  such  parameters  as  tracking  aperture  size, 
target  motion,  and  diversity  techniques.  The  entry  of  the  amplitude 
and  phase  data  for  the  scattering  complexes  was  implemented  through 
a  separate  subroutine  to  facilitate  use  of  the  programs  with  measured 
static  range  patterns,  should  they  become  available.  Detailed  flow 
diagrams  and  listings  of  the  various  programs  and  subroutines  which 
evolved  are  presented  in  Appendix  B. 

A,  Program  Philosophy 

Fig.  5-1  is  a  summary  flow  diagram  which  illustrates  the  basic 

computer  implementation.  The  parameters  which  are  read  into  the 

program  at  execution  time  include  radar  wavelength,  aperture  diameter, 

coordinates  and  identifying  codes  for  the  various  scattering  complexes, 

target  motion  data,  and  the  sampling  time  increment.  The  space 

trajectory  over  which  glint  or  other  field  parameters  are  to  be 
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determined  is  established  by  a  subroutine  called  TRAJ.  Any  desired 
trajectory  is  thus  easily  incorporated  simply  by  preparing  a  subroutine 
which  describes  the  desired  series  of  observation  points,  in  spherical 
coordinates,  as  a  function  of  time.  The  field  parameters  are  then 
computed  at  each  of  these  sampling  points. 

The  first  operation  in  the  computation  of  the  field  parameters 
is  the  appropriate  modification  of  the  coordinates  obtained  from  TRAJ 
to  account  for  any  change  of  aspect  angle  due  to  target  motion.  The 
modified  angular  coordinates  of  the  observation  point  are  then 
transformed  into  equivalent  coordinates  with  respect  to  each  scattering 
complex,  and  these  values  used  to  select  the  appropriate  amplitude 
and  phase  values  from  the  scattering  complex  data.  This  selection  is 
provided  by  a  subroutine  named  AMPPHI,  which  may  be  written  to  provide 
either  values  for  a  desired  group  of  theoretical  scatterers,  or  values 
from  actual  recorded  data  from  scattering  complexes.  The  resulting 
phase  values  are  then  adjusted  to  include  any  Doppler  effects  due  to 
target  motion. 

Haying  obtained  the  appropriate  amplitude  and  phase  values  for 
each  complex,  the  next  computation  involves  the  determination  of  the 
distances  from  the  observation  point,  P,  to  each  complex,  and  the 
distances  from  each  complex  to  the  extremes  of  the  aperture.  These 
propagation  paths  are  then  used  to  find  the  field  contributions  of 
each  complex  at  each  extreme  of  the  aperture,  and  the  contributions 
summed  to  obtain  the.  total  field  values. 

The  glint  in  each  tracking  plane  is  determined  from  the  total 


field  phase  at  the  aperture  extremes  as  indicated  in  equation  (4-30) . 

The  time,  glint  error  in  each  plane,  total  returned  amplitude  and 
phase  at  P,  and  the  amplitude  contribution  at  P  due  to  each  complex 
are  written  into  magnetic  disk  storage  for  future  processing  and 
analysis.  The  program  then  transfers  control  to  repeat  the  computations 
at  the  next  sample  point  along  the  trajectory. 

B.  Synthesis  of  Scattering  Complex  Patterns 
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summation  of  the  contributions  from  many  different  spheres.  Because 
many  different  data  points  were  required  to  describe  adequately  each 
pattern,  it  became  apparent  that  the  computation  time  for  this  synthesis 
technique  was  quite  large.  For  this  reason,  a  more  direct  analytical 
approach  was  needed. 

Another  pattern-generation  technique  which  could  be  used  is  to 
describe  the  desired  function  of  aspect  angle  by  an  algebraic  equation. 
An  equation  of  the  form 

f  C4>)  =  K  cos(n4>),  (5-1) 

where  K  is  a  constant  and  n  is  an  integer,  generates  a  polar  lobe 
pattern.  As  n  becomes  large,  the  number  of  lobes  also  increases.  The 
product  of  two  functions  of  this  form  can  be  used  to  generate  a  large 
variety  of  polar  patterns.  Expressed  as  a  logarithmic  (power)  amp¬ 
litude  function,  such  a  product  is 

2 

A  =  10  log  {K  cos(n  d>)  cos(n7<j>)]  dB.  (5-2) 

10  1  ^ 

Here  K  is  some  maximum  amplitude,  is  some  small  integer,  and  v-2  is 
a  large  integer.  Thus  the  latter  cosine  function  provides  the  desired 
degree  of  fine-grain  lobe  structure,  while  the  former  results  in  a 
slowly-varying  amplitude  envelope.  Fig.  5-2  is  a  plot  of  equation 
(5-2)  with  K  =  10,  n^  =  1  and  ^  =  90.  A  phase  pattern  can  be  generated 
in  a  similar  manner  by  use  of  a  function  of  the  form 


(5 


a  -  180  cos(n^<j>)  deg 

where  the  integer  is  chosen  to  provide  the  desired  rapidity  of 
phase  change  with  aspect  angle. 

Because  of  the  relative  simplicity  of  this  second  technique  of 
pattern  synthesis,  it  was  used  in  most  of  the  computer  studies.  A 
further  simplification  was  to  restrict  trajectories  used  with  a 
synthesized  target  to  elevation  angles  near  0  =  90° ,  near  the 

azimuth  plane..  Table  5-1  summarizes  the  amplitude  and  phase 
generating  functions  for  four  synthesized  scattering  complexes  which 
were  used  in  the  studies.  The  additional  sine  functions  involving 
0  were  included  in  the  amplitude  equations  to  provide  a  degree  of 
variation  in  the  elevation  plane.  Table  5-2  lists  the  spherical 
coordinates  of  each  of  these  four  complexes,  relative  to  a  primary 
origin,  which  were  used  for  many  of  the  computer  data  runs.  This 
arrangement  was  assigned  the  designation  configuration  No.  10. 

G.  Illustrations  of  Phase  Front  Behavior 
To  provide  insight  into  the  nature  cf  phase  front  distortion, 
it  is  instructive  to  examine  the  phase  front  patterns  for  several 
target  configurations.  These  can  be  obtained  by  describing  a 
trajectory  which  varies  over  the  desired  region  of  aspect  angles 
while  maintaining  a  constant  range,  converting  the  computed  total 
field  phase  at  each  sample  point  to  the  equivalent  phase  front  by  use 
of  ,  and  plotting  the  result. 
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TABLE  5-1 

AMPLITUDE  AND  PHASE  GENERATING  FUNCTIONS  FOR  SYNTHESIZED 
SCATTERING  COMPLEXES  20,  21,  22,  23 


COMPLEX 


GENERATING  FUNCTIONS 


L, 


20 


A  =  10  log j^q  [sin  (0)  sin  ($)  cos  (360  <j>)]^ 
a  =  180  cos  (720  <J>) 


21 


22 

& 

23 


A  =  10  log^  [0.9  sin  (6)  sin  (<j>)  cos  <  //  <}>)  ] ^ 

a  =  180  cos  (360  (j>) 

2 

A  =  10  log^Q  [sin  (0)  cos  (<j>)  cos  (360  4>) ] 

(—it / 2  <  d  <  tt/2) 

A  =  0  (<t>  >  n/2) 


L 


} 


a  =  180  cos  (720  4>) 


6.1 


Complex 


TABLE  5-2 


SPHERICAL  COORDINATES  FOR  SCATTERING  COMPLEXES 
IN  TARGET  CONFIGURATION  NO.  10 


d  (m) 


<j>'(deg)  6"(deg) 
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Fig.  5-3  is  a  plot  of  two  different  phase  front  patterns  for 
two  isotropic  reflectors.  Note  that  the  ordinate  is  plotted  full- 
scale  in  units  of  inches.  The  abcissa  is  the  azimuth  angle,  <j>,  as 
measured  from  the  line  joining  the  reflectors.  In  the  upper  plot, 
the  amplitude  ratio  between  the  reflectors  is  0.9,  while  the  lower 
pattern  was  obtained  tor  a  ratio  of  1.5.  In  both  cases  the  spacing 
was  one  wavelength  at  10  GHz,  and  the  range  was  1  Km.  The  apparent 
direction  to  the  midpoint  of  the  reflectors  from  any  given  observation 
point  is  given  by  the  normal  to  the  phase  front  at  that  point.  From 
the  figure,  it  is  seen  that  the  maximum  slope  of  the  phase  front  is 
much  larger  for  the  amplitude  ratio  nearest  to  unity  (A~0„ 9) ,  which 
results  in  a  larger  glint  error.  Also,  the  magnitude  of  the  ratio  as 
it  compares  to  unity  governs  the  algebraic  sign  of  the  error,  jL.e,. 
to  which  side  of  the  mid-point  of  the  reflectors  the  apparent  target 
center  lies.  Both  of  these  observations  are  compatible  with  the 
results  of  Meade,  as  depicted  in  Fig.  2-3. 

A  full-sea1 e  phase  front  plot  for  the  synthesized  complex 
target  which  is  described  in  Tables  5-1  and  5-2  is  presented  in  Fig. 
5-4.  The  azimuth  angle  is  measured  from  the  positive  x-axis  (4>*=0) , 
and  the  observations  are  made  in  the  0=90°  plane  at  a  range  of  1  Km. 
The  radar  frequency  is  10  GHz.  Note  that  the  larger  size  of  the 
target  configuration  produces  a  muc  .•  re  rapid  variation  of  the 
phase  front  with  aspect  angle  than  that  shown  in  Fig.  5-3. 


PHASE  FRONTS  FOR  TWO  REFLECTORS  SPACED  ONE  WAVELENGTH  APART 
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Comparison  with  results  of  Mead 


A  glint  error  determination  was  made  for  the  two-reflector  target 
used  in  the  example  problem  in  Chapter  2.  The  trajectory  consisted 
of  a  single  sample  point  at  a  range  of  1000  m  along  the  perpendicular 
to  a  line  joining  the  two  reflectors.  An  AMPPHI  subroutine  was  written 
which  incremented  the  phase  for  one  of  the  reflectors  in  5  degree 
steps  and  also  varied  the  amplitude  ratio  >  A,  between  the  reflectors. 
The  reflectors  were  spaced  one  wavelength  at  the  10  GHz  radar 
frequency,  and  the  glint  error  for  each  combination  of  relative  phase 
and  amplitude  was  normalized  to  this  spacing.  Fig.  5-5  is  a  family  of 
curves  depicting  the  results  for  various  values  of  amplitude  ratio. 

The  figure  compares  quite  favorably  with  Fig.  2-3;  in  fact,  the  computer 
generated  plots  appear  to  be  overlays. 

Range  dependence  of  glint 

If  linear  glint  is  independent  of  range,  an  angular  boresight¬ 
tracking-error  value  from  a  radar  measurement  could  be  converted  to  a 
measurement  independent  of  range  by  multiplying  the  value  by  the 
range  to  the  target.  Such  a  range  dependence  for  two  isotropic 
reflectors,  and  R^,  can  be  illustrated  by  consideration  of  Fig.  5-6. 
The  origin  is  designated  by  0,  and  the  apparent  target  center  by  C. 

A  aeries  of  observation  points  lie  along  the  radial  line  OP'.  The 
arcs  centered  on  0  and  C  represent  the  phase  fronts  which  would 
result  from  a  single  reflector  at  0  or  C,  respectively.  Note  that 


this  construction  yields  a  phase-front -difference  between  the  pairs 


of  arcs  which  is  constant  with  range  along  the  radial  lines  OA  and  OB. 
This  result  is  compatible  with  the  concept  of  phase  retardation  of  a 
propagating  wave.  Note  also  that  the  angular  glint  decreases  with 
range.  From  a  physical  viewpoint,  this  means  that  at  a  sufficiently 
large  range,  the  target  configuration  appears  as  a  single  point  which 
produces  amplitude  scintillation  but  induces  no  appreciable  angular 
glint.  Fig.  5-7  shows  the  linear  and  angular  glint  dependence  for 
two  isotropic  reflectors  spaced  two  meters  apart  at  10  GHz. 

For  the  case  of  a  target  configuration  more  complex  than  that 
just  described,  it  appears  that  the  range  dependence  can  no  longer  be 
easily  described.  The  different  combinations  of  the  various  target 
elements  produce  lobing  patterns  which  have  effective  centers  away 
from  the  origin.  Thus  changing  range  changes  the  relative  contributions 
from  the  various  scattering  centers,  with  a  resulting  change  in  the 
linear  glint  error.  To  illustrate  this  effect,  the  two  isotropic 
reflectors  used  in  generating  Fig.  5-7  were  replaced  by  two  identical 
nonisotropic  reflectors  having  amplitude  patterns  similar  to  Fig.  5-2. 
The  resulting  glint  behavior  with  range  is  shown  in  Fig.  5-8.  Note 
that  the  angular  glint  continues  to  show  a  decreasing  trend  with 
increasing  range,  in  accordance  with  the  physical  reasoning  discussed 
earlier. 

Analysis  of  glint  waveforms  for  a  synthesized  target 

A  plot  of  linear  glint  error  vs  time  for  configuration  No.  10 
is  shown  in  Fig.  5-9.  The  component  of  error  in  the  transverse 
tracking  plane  was  computed  and  plotted  at  0.025  second  increments 
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along  the  trajectory  listed  in  Appendix  B  as  trajectory  No.  11.  This 
trajectory  provides  a  series  of  sample  points  having  a  maximum  range 
of  4.65  Km,  a  minimum  range  of  900  m,  and  aspect  angles  which  vary 
about  nominal  values  of  8=90°,  i^S0.  The  projected  target  extent 
from  this  aspect  is  approximately  20  meters.  The  root-mean-squared 
(RMS)  value  of  this  waveform  is  11.9  meters,  as  measured  from  the 
origin. 

Fig.  5-10  presents  a  probability  distribution  for  the  glint 
waveform  of  Fig.  5-9.  This  plot  was  generated  by  a  program  which 
divided  the  total  range  of  glint  error  shown  on  the  abscissa  into 
one  meter  intervals  and  scanned  the  1001  data  points  of  Fig.  5-9 
to  determine  the  number  of  values  falling  within  each  interval.  Each 
of  these  numbers  was  then  normalized  to  the  total  number  of  data 
points  and  plotted  vs^  the  abscissa  value  corresponding  to  the  midpoint 
of  the  appropriate  interval. 

The  spectral  density  for  the  waveform  of  Fig.  5-9  is  shown  in 
Fig.  5-11,  No  attempt  was  made  to  calibrate  the  spectral  determinations 
in  terms  of  meters/>/Sz’,  therefore  the  ordinate  values  can  be  used  only 
as  a  relative  indication  of  the  frequency  content  of  the  glint  waveform. 
The  published  experimental  curves  of  Howard  and  others  [2,3,8,24,29] 
indicate  little  spectral  energy  above  about  15  Hz.  This  value  was 
increased  to  20  Hz  as  a  maximum  expected  frequency  component,  and 
Shannon’s  sampling  theorem  used  to  establish  the  sampling  increment 
of  0.025  seconds.  The  25  seconds  of  available  data  estaolished  the 
lowest  resolvable  frequency  at  0.04  Hz,  The  plotted  values  are  the 
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results  of  a  Fourier  analysis  at  0.04  Hz  and  every  multiple  of  0.04 
Hz,  up  to  a  maximum  of  7  Hz.  Assuming  the  upper  cutoff  of  a  typical 
tracking  servo  bandpass  to  be  1  Hz,  most  of  the  energy  in  the  wave¬ 
form  would  enter  the  system  as  tracking  noise. 

As  an  illustration  of  the  change  in  the  glint  waveform  for  a 
different  aspect  angle.  Fig.  5-12  shows  the  glint  error  vs_  time 
along  trajectory  No.  15.  This  trajectory  differs  from  No.  11  only  in 
the  choice  of  $=48°  as  the  nominal  azimuth  aspect  angle.  Fig.  5-13 
and  Fig.  5-14  show  the  probability  distribution  and  spectral  density, 
respectively,  for  Fig.  5-12.  Note  that  although  the  RMS  and  most 
probable  values  of  apparent  target  position  have  changed,  most  of 
the  larger  spectral  components  remain  below  1  Hz. 

Relationship  between  glint  and  field  amplitude 

It  may  be  recalled  from  the  discussion  in  Chapter  III  that  several 
investigators  indicated  a  negative  correlation  between  regions  of  large 
glint  error  and  total  field  amplitude.  This  hypothesis  was  tested 
using  synthesized  configuration  No.  10  along  two  different  trajectories 
and  at  several  different  radar  frequencies.  For  each  case,  the  linear 
correlation  between  the  total  field  amplitude  and  the  magnitude  of  the 
linear  glint  error  was  computed.  The  results  presented  in  Table  5-3 
show  the  correlation  coefficients  to  be  negative  and  quite  large  -  - 
significant  indeed  for  correlation  between  two  different  types  of 
field  parameters.  This  strong  negative  correlation  is  illustrated 
graphically  by  Fig.  5-15,  in  which  the  glint  error  waveform  of  Fig.  5-9 


is  superimposed  with  a  relative  plot  of  the  corresponding  total  field 
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amplitude.  Note  that  the  regions  of  large  positive  and  negative  error 


coincide  exactly  with  the  regions  of  small  amplitude.  The  linear 
correlation  between  amplitude  and  magnitude  of  glint  error  for  this 
figure  is  -0.605. 

Effect  of  target  yaw 

In  order  to  illustrate  the  effect  of  constant  angular  motion  of 
the  synthesized  target  on  glint  behavior,  subroutines  were  prepared 
to  include  the  appropriate  coordinate  and  phase  modifications  for  the 
case  of  target  yaw  with  an  observation  point  in  the  azimuth  plane. 

A  plot  of  glint  error  vs_  time  for  a  yaw  rate  of  0.0005  rad/sec  is 
presented  in  Fig.  5-16.  With  the  exception  of  the  non-zero  yaw  rate, 
all  parameters  are  identical  to  those  producing  Fig.  5-9.  The 
linear  correlation  coeffecient  between  the  sets  of  error  values  shown 
in  these  two  figures  is  -0.14,  indicating  an  effective  decorrelation 
of  the  glint  dat'..  Probability  distribution  and  spectral  density  plots 
for  the  waveform  of  Fig,  5-16  appear  in  Fig.  5-17  and  Fig.  5-18, 
respectively.  From  these  figures,  it  appears  that  although  the  waveform 
with  target  yaw  is  different,  the  characteristics  of  the  waveform 
remain  quite  similar.  The  RMS  error  value  is  essentially  unchanged 
while  the  probability  distribution  still  shews  significant  values  to 
+  25  meters.  The  spectral  content  has  been  shifted  slightly  towards 
the  lower  frequencies. 

Influence  of  scattering  range  angular  data  increment 

An  ultimate  limitation  in  the  usefulness  of  the  model  may  be  the 
accuracy  with  which  scattering  range  pattern  recorders  can  measure 
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the  angular  setting  of  the  target  positioner.  The  required  accuracy, 
of  course,  will  depend  upon  the  fine-grain  lobe  structure  of  the 
scattering  complexes.  To  determine  the  effect  of  incrementing  the 
aspect  angle  for  a  synthesized  target  configuration,  an  AMPPHI 
subroutine  was  written  which  not  only  implemented  the  generating 
functions  listed  in  Table  5-1,  but  also  allowed  for  rounding  the 
values  of  the  aspect  angles  to  be  used  in  generating  the  amplitude 
and  phase  data.  That  is,  the  aspect  angles  of  each  observation 
point  in  each  secondary  coordinate  system  were  rounded  to  the 
nearest  "Ru  before  being  used  to  obtain  the  appropriate  amplitude 
and  phase  data.  The  result  of  a  linear  correlation  between  the 
error  data  of  Fig.  5-9  and  an  identical  run  with  R=0.25  deg  was  a 
correlation  coefficient  of  only  0.071.  This  seems  to  indicate  that, 
at  least  for  complexes  with  lobe  structures  similar  to  those  of  the 
synthesized  scattering  complexes,  a  very  accurate  pattern  range 
positioner  would  be  required. 


E,  Diversity  Effects 

It  was  implied  at  the  beginning  of  Chapter  IV  that  the  intended 
utility  of  the  glint  model  developed  in  that  chapter  was  primarily 
one  of  providing  a  means  of  testing  and  evaluating  various  implements 
tion  of  glint  compensation  schemes.  As  an  illustration  of  this  use, 
this  last  section  will  demonstrate  the  effectiveness  of  two  simple 
data  processing  schemes  for  dual  frequency  and  dual  tracker  diversity 
techniques.  The  target  used  is  that  which  has  been  designated  as 


synthesized  configuration  No.  10. 


Frequency  diversity 

Table  5-4  presents  typical  correlation  coefficients  between  sets 
of  error  data  obtained  over  trajectory  No.  11  with  frequency  as  a 
parameter.  A  correlation  of  0.5  or  less  was  chosen  as  a  reasonable 
condition  for  independence  of  the  data.  Ray  [32]  has  used  a  similar 
value.  Thus,  judging  from  the  values  in  Table  5-4,  it  would  seem 
advantageous  to  investigate  the  use  of  the  glint  data  obtained  at 
two  or  more  of  these  frequencies  to  obtain  more  accurate  tracking 
information.  The  magnitude  of  RMS  error  before  and  after  diversity 
processing  will  be  used  as  a  measure  of  accuracy  improvement  (or 
deterioration) . 

The  first  technique  implemented  was  to  simply  average  the  two 
values  obtained  with  the  two  frequencies  at  each  sample  point  along 
the  trajectory.  The  results  of  this  operation  are  presented  in 
Table  5-5,  where  the  RMS  values  listed  are  in  meters,  and  the  values 
along  the  diagonal  are  the  RMS  values  at  the  corresponding  single 
frequencies.  Note  that  in  most  cases,  the  result  of  averaging  at 
two  frequencies  is  a  value  less  than  the  value  at  either  single 
frequency. 

Another  technique  of  data  processing  which  was  implemented  was 
based  upon  the  negative  correlation  between  glint  error  and  total 
field  amplitude  demonstrated  earlier  in  this  chapter.  The  computer 
program  selected  from  the  pair  of  glint  values  at  each  sample  point 
that  value  which  corresponded  to  the  frequency  having  the  largest 
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field  amplitude.  The  resulting  RMS  error  values  are  presented  in 
Table  5-6,  where  the  values  on  the  diagonal  are  again  the  single¬ 
frequency  RMS  errors.  This  processing  technique  appears  to  provide 
improved  tracking  over  simple  averaging. 

To  illustrate  the  nature  of  the  improvement  provided  by 
diversity,  it  is  instructive  to  analyze  a  resulting  glint  waveform 
and  compare  it  with  the  corresponding  waveforms  at  the  single 
frequencies.  The  frequency  combination  in  Table  5-6  resulting  in 
the  lowest  RMS  value  is  10  -  9.9  GHz.  The  waveform,  probability 
distribution,  and  spectral  density  plots  for  f  =  10  GHz  have  been 
presented  in  Figures  5-9,  5-10  and  5-11.  The  corresponding  plots 
for  f  =  9.9  GHz  are  shown  in  Figures  5-19,  5-20  and  5-21.  The 
graphs  resulting  from  using  the  error-selection  technique  at  these 
two  frequencies  appear  as  Figures  5-22,  5-23  and  5-24.  Note  that 
the  large  values  of  error  are  eliminated,  with  a  corresponding 
narrowing  of  the  probability  distribution  function  about  zero.  Also, 
the  spectral  lines  seem  to  be  more  concentrated  near  zero  frequency. 

As  a  test  of  a  slightly  more  sophisticated  diversity  processing 
scheme,  the  error-selection  technique  was  implemented  on  the  basis 
of  the  amplitudes  at  all  six  of  the  radar  frequencies  previously 
used.  This  resulted  in  a  RMS  value  of  only  4.6  meters. 

Space  diversity 

Perhaps  the  most  simple  form  of  space  diversity  is  the  phase  front 
averaging  effect  provided  by  a  finite-sized  antenna  aperture.  All  of 
the  data  presented  thus  far  in  this  chapter  has  been  generated  with  an 
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aperture  diameter  of  50  cm.  To  gain  insight  concerning  the  degree  of 
aperture  averaging  which  might  be  obtainable,  two  runs  with 
configuration  No.  10  along  trajectory  No.  11  and  at  10  GHz  were  made, 
one  with  a  diameter  of  10  cm  and  the  other  using  a  diameter  of  100 
cm.  The  linear  correlation  coefficient  between  the  glint  errors  for 
these  two  cases  was  0.96.  This  tends  to  indicate  that  even  an  order 
of  magnitude  of  change  in  aperture  diameter  yields  little  decorrelation 
of  the  data. 

The  more  usual  form  of  space  diversity  involves  the  use  of  two 
tracking  antennas  spaced  some  distance  apart.  While  this  technique  is 
not  practical  in  a  homing  missile  application,  several  runs  were  made 
for  configuration  No.  10  at  a  frequency  of  10  GHz  and  using  a  dual 
tracker  to  gain  insight  concerning  the  degree  of  decorrelation  achieved 
by  spaced  antennas.  Each  dual  tracker  run  was  actually  implemented 
by  making  two  data  runs  of  the  main  glint  program.  In  each  case,  the 
TRAJ  subroutine  generated  the  normal  trajectory  No.  11,  but  for  each 
observation  point,  the  value  of  the  azimuth  coordinate  was  increased 
or  decreased  by  an  amount  corresponding  to  a  fixed  linear  distance 
at  the  value  of  range  for  the  observation  point.  A  linear  correlation 
was  then  performed  on  the  pairs  of  glint  error  values  obtained  for 
the  two  data  runs,  and  RMS  computations  made  for  each  data  set  and 
also  on  the  result  of  using  the  error-selection  technique  based  upon 
the  received  amplitude  at  each  tracker.  The  results  of  these 
computations  are  presented  in  Table  5-7,  where  FMS^  and  RMS^  are  the 
RMS  error  values  for  the  individual  trackers,  and  KMD,.  is  the 


corresponding  value  using  error-selection  diversity.  It  appears 
from  the  data  that  even  rather  small  spacings  can  be  quite  effective 
in  reducing  glint  error. 

There  is  another  form  of  space  diversity  which  is  more  practical 
in  a  homing  missile  application.  This  involves  the  use  of  bistatic 
tracking  with  two  or  more  target  illuminating  antennas.  Because 
bistatic  target  illumination  is  often  ground  situated,  the  use  of 
multiple  antennas  at  larger  spacings  becomes  practical.  The  basic 
program  limitation  to  the  monostatic  configuration  prevented 
investigation  of  this  technique,  but  of  course  the  model  as  developed 
could  be  adapted  for  bistatic  studies. 


V.  CONCLUSIONS 


In  summary,  a  review  of  prior  glint  investigations  reveals 
no  analytical  model  which  is  adequate  in  describing  glint  for  a 
complex  target,  and  yet  simple  enough  to  permit  implementation  for 
such  a  realistic  scattering  object.  The  technique  developed  here 
uses  the  concept  of  coherent  summation  of  the  field  contributions 
from  various  nonisotropic  scattering  complexes.  The  glint  error  is 
determined  by  use  of  a  simple  geometric  construction  together  with 
values  of  total  field  phase  in  the  region  of  the  tracking  aperture. 

The  model  was  implemented  on  the  digital  computer,  and  the 
results  obtained  for  theoretical  configurations  of  isotropic 
scatterers  were  analyzed  and  compared  with  known  results  from  other 
investigations.  In  lieu  of  measured  data  from  actual  scattering 
complexes,  nonisotropic  amplitude  and  phase  patterns  were  synthesized 
to  illustrate  further  the  modelling  technique.  The  glint  waveforms 
for  theoretical  target  configurations  of  such  synthesized  complexes, 
as  obtained  over  several  different  approach  trajectories  and  at 
several  different  radar  frequencies,  were  analyzed  for  RMS  value, 
probability  distribution,  and  spectral  content. 

Finally,  the  effects  of  using  frequency  and  space  diversity  with 
a  synthesized  target  were  studied.  The.  results  indicate  that  a 
considerable  tracking  improvement  can  be  achieved  by  using  the  easily 
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implemented  data  processing  scheme  of  choosing  the  tracking-error 
signal  which  corresponds  to  the  diversity  return  having  the 
largest  signal  amplitude. 

The  data  obtained  and  the  results  discussed  indicate  that  thii 
study  could  be  profitably  pursued  both  to  gain  valuable  insight, 
into  the  nature  of  glint,  and  as  a  tool  to  devise  and  evaluate 
diversity  compensation  schemes.  It  seems  that  this  end  would 
adequately  justify  the  time  and  expense  of  obtaining  the  required 
amplitude  and  phase  patterns  for  actual  scattering  complexes. 


REFERENCES 


1.  Locke,  A.  S.,  Guidance.  Princeton,  N.  J.,  D.  Van  Nostrand,  Co., 
Inc.  (1955). 

2.  Dunn,  J.  H.,  Howard,  D.  D.,  and  King,  A.  M.,  "Phenomena  of 
Scintillation  Noise  in  Radar-Tracking  Systems,"  Proe.  IRE,  47, 

855  (1959). 

3.  Meade,  J.  L.,  Hastings,  A.  E.,  and  Gerwin,  H.  L.,  Noise  in  Tracking 
Radars.  Naval  Research  Laboratory,  Report  3759,  AD  90  876  (1950). 

4.  BrockneF,  C.  E.,  "Angular  Jitter  in  Conventional  Conical  Scanning 
Automatic  Tracking  Radar  Systems,"  Proc.  IRE,  39,  51  (1951). 

5.  Delano,  R.  H.,  "A  Theory  of  Target  Glint  or  Angular  Scintillation 
in  Radar  Tracking,"  Proc.  IRE,  41,  1778  (1953). 

6.  Kerr,  D.  E.,  Propagation  of  Short  Radio  Waves.  New  York,  N.  Y., 
McGraw-Hill  Book  Co.  (1951). 

7.  Muchmore,  R.  B.,  "Aircraft  Scintillation  Spectra,"  IRE  Trans.  AP, 

8,  201  (1960). 

8.  Muchmore,  R.  B.,  Review  of  Scintillation  Measurements.  Hughes 
Aircraft  Co.  TM  272,  AD4849  (1952). 

9.  Watson,  R.  B.,  and  Feagin,  R.  B. ,  Analysis  of  Angular  Errors  of 
a  Track  Radar  System  for  N-Point  Radar  Target.  Military  Physics 
Research  Lab.,  University  of  Texas  (1956). 

10.  Feagin,  R.  B. ,  Watson,  R.  B.,  and  Norwood,  J.  M. ,  Detailed  Data 
Reduction  Procedure  for  Airborne  N-Point  Target  Backscattering. 
Military  Physics  Research  Lab.,  University  of  Texas  (1955). 

11.  Weimer,  F.  C. ,  and  Peters,  Leon,  Jr.,  Study  of  Pointing  Errors  in 
Conically  Scanning  and  Monopulse  Tracking  Radars  for  Multi-Point 
Targets.  Antenna  Lab.,  Ohio  State  University,  601-12,  AD  123  476 
(1956). 

12.  Peters,  Leon,  Jr.,  Accuracy  of  Tracking  Radar  Systems.  Antenna 
Lab.,  Ohio  State  University,  601-29.  AD  200  027  (1957). 

13.  Lewis,  D.  E.,  On  Tracking  Radar  Systems  and  the  Target  Phase 
Center.  Antenna  Lab,,  Ohio  State  University,  777-4,  AD  211  187 
(1958). 


103 


104 

Howard,  D.  D.,  Radar  Target  Angular  Scintillation  tn  Tracking  and 
Guidance  Systems  Based  on  Echo  Signal  Phase  Front  Distortion. 

Naval  Research  Laboratory,  Report  5430  (1959). 

Goldstein,  M.  H.  Jr.,  and  Brindey,  W.  E.,  Glint  Noise  and  Signal 
Simulation  for  the  Meteor  Seeker.  MIT  Research  Memo  RM  6978-3, 

AD  58  119  (1952). 

Delano,  R.  H.,  and  Pfeffer,  I.,  "The  Effect  of  AGC  on  Radar 
Tracking  Noise,"  Proc.  IRE,  44,  801  (1956). 

Dunn,  J.  H.,  and  Howard,  D.  D. ,  "The  Effects  of  Automatic  Gain 
Control  Performance  on  the  Tracking  Accuracy  of  Monopulse  Radar 
Systems,"  Proc.  IRE,  47,  430  (1959). 

Peters,  Leon,  Jr.,  and  Weimer,  F.  C.,  "Concerning  the  Assumption 
of  Random  Distribution  of  Scatterers  as  a  Model  of  an  Aircraft 
for  Tracking  Radars,"  IRE  Trans.  AP,  9,  110  (1961). 

Muchmore,  R.  B.,  "Reply  to  Comments  by  Leon  Peters,  Jr.  and 
F.  C.  Weimer,"  IRE  Trans.  AP,  9,  112  (1961). 

Delano,  R.  H. ,  "Further  Reply  to  Comments  by  Leon  Peters,  Jr. 
and  F.  C.  Weimer,"  IRE  Trans.  AP,  9,  227  (1961). 

Peters,  Leon,  Jr.,  and  Weimer,  F.  C.,  "Reply  to  Comments  by 
R.  H.  Delano,"  IRE  Trans.  AP,  9,  228  (1961). 

Peters,  Leon,  Jr.,  and  Weimer,  F.  C.,  "Reply  to  Comments  by 
R.  B.  Muchmore,"  IRE  Trans.  AP,  9,  229  (1961). 

Siegel,  K.  M. ,  "Comparison  Between  Theoretical  and  Experimental 
Radar  Cross  Sections  of  Aircraft,"  IRE  Trans.  AP,  9,  228  (1961). 

Povejsil,  D.  J.,  Raven,  R.  S.,  and  Waterman,  P. ,  Airborne  Radar . 
Princeton,  N.  J. ,  D.  Van  Nostrand,  Co.,  Inc.  (1961). 

Moll,  J.  W.,  "On  the  Radar  Echo  from  Aircraft,"  IEEE  Trans .  on 
Aerospace  and  Elec,  Systems,  574  (1967) . 

Ostrovityanov,  R.  B.,  "Angular  Noise,"  Radiotekhnika  i  Elektronika, 
4,  307  (1966) . 

Gubonin,  N.  S.,  "Fluctuation  of  the  Phase  Front  of  the  Wave 
Reflected  from  a  Complex  Target,"  Radiotekhnika  i  Elektronika, 

10,  717  (1965). 

Coleman,  J.  T.,  and  Someroski,  F.  W. ,  A  Survey  and  Analysis  of 
Radar  Glint  As  Applied  to  Air  Defense  Missile  System  Design. 
Battelle  Memorial  Institute,  RSIC-699,  AD  385  494  (1967). 


105 

29.  Dunn,  J.  H.,  and  Howard,  D,  D.,  "Radar  Target  Amplitude,  Angle, 
and  Doppler  Scintillation  from  Analysis  of  the  Echo  Signal 
Propagating  in  Space,"  IEEE  Trans.  MTT,  16,  715  (1968). 

30.  Birkemeir,  W.  P.,  and  Wallace,  N.  D.,  "Radar  Tracking  Accuracy 
Improvement  by  Means  of  Pulse-to-Pulse  Frequency  Modulation," 

AIEE  Trans.  Comm,  Elec,,  81,  571  (1963) . 

31.  Gaheen,  A.  F.,  McDonough,  J.  T.,  and  Tice,  D.  P.,  Frequency 
Diversity  Radar  Study:  Vol  I.  Westinghouse  Defense  and  Space 
Center ,  AFAL-Tk-66-128,  AD  480  358  (1966). 

32.  Ray,  H.  K. ,  "Improving  Radar  Range  and  Angle  Detection  with 
Frequency  Agility,"  Microwave  J.,  9,  63,  (1966). 

33.  Senior,  T.  B.  A.,  The  Elevation  Errors  Involved  in  the  Radar 
Location  of  a  Low  Flying  Target.  Royal  Radar  Establishment, 

Memo  1220  (1956). 

34.  Shaw,  A.  H.,  An  Expression  for  the  Ground  Reflection  Elevation 

Error  of  an  Auto-Tracking  Radar  at  Elevation  Angles  Comparable 
with  the  Beamwidth.  Royal  Radar  Establishment,  Memo  1457  (1961). 

35.  Graf,  E.  R.,  Weathers,  G.  D.,  Sim>,  R.  J.,  and  Fincher,  C.  F., 

A  Study  of  Radar  Tracking  Error  Resulting  from  Multipath:  Analysis 
and  Experimental  Results.  Auburn  University  Engineering  Experiment 
Station,  Tech.  Rept,  No.  5,  AMC- 14525  (1967). 

36.  RAT  SCAT  Radar  Target  Scatter- Site .  Holloman  AFB,  New  Mexico, 
MDC-TR-65-26  (1965) 

37.  Caster,  J.  M.,  and  Walier,  A. -A.,  Jr.,  Significance  of  Target 
Characteristics  upon  Precision  of  Monopulse  Radar.  Syracuse 
University  Research  Jorp.,  DSL  R-94,  AD  451  663  (1964). 

38.  Allen,  P.  J.,  and  Olin,  I.  D. ,  A  Four- Component  Polarization 
Resolver  for  Microwaves.  Naval  Research  Lab.,  Report  1086, 

AD  646  189  (1960). 

39.  Olin,  I.  D.,  A  System  for  the  Dynamic  Measurement  of  Radar  Return 
Polarization.  Naval  Research  Lab.,  Report  1082,  AD  646  188  (1960). 

40.  Olin,  I.  D.,  and  Queen,  F.  D.,  Measurements  Using  a  Polarization 
Instrumentation  Radar  on  Navigational  Buoys.  Naval  Research  Lab., 
Report  5701,  AD  268  727  (1961). 

41.  Berkowitz,  R.  S.,  Modern  Radar.  New  York,  N.  Y.,  John  Wiley  and 
Sons,  Inc.  (1965). 


106 

42.  Graves,  C.  D.,  Radar  Polarization  Power-Scattering  Matrix. 
University  of  Michigan,  2144-45-T,  AD  103  803  (1955). 

43.  Ivey,  H.  D.,  and  Long,  M.  W.,  Polarization  Characteristics  of 
Radar  Targets.  Georgia  Institute  of  Technology,  A-158-1, 

AD  649  58  (1954) . 


rr 


107 


APPENDIX  A 

COORDINATE  AND  PHASE  MODIFICATIONS  FOR  UNIFORM 
ROLL,  PITCH,  AND  YAW  RATES 


The  purpose  of  this  appendix  is  to  develop  the  appropriate  coordi¬ 
nate  and  phase-term  modifications  for  the  case  of  uniform  target  rretion 
about  the  roll,  pitch,  and  yaw  (x,y,z,)  axes.  Uniform  yaw  will  be 
examined  first,  followed  by  a  discussion  of  the  required  adaptations 
to  include  roll  and  pitch. 

Because  yaw  involves  rotation  of  the  target  about  the  z-axis,  the 
change  in  aspect  angle  can  be  included  simply  by  modifying  the  azimuth 
(<{>)  coordinate  of  the  aspect  point  P.  If  the  elapsed  time  since  some 
reference  time  tQ  is  T  seconds,  and  the  yaw  rate  is  toy  rad/sec,  the 
modified  azimuth  coordinate  of  P  is 

4»  *  “  <}>  +  WyT  (A-l) 

where  is  taken  to  be  positive  in  the  direction  of  increasing  <j>. 

In  addition  to  the  effective  change  in  position  of  P,  the  yawing 
motion  also  causes  a  Doppler  shift  in  the  frequency  of  the  radar  return 
from  each  complex  C^.  This  effect  can  be  included  by  making  an 
appropriate  modification  in  the  phase  of  the  signal  contribution  from 
e  <  i  complex  before  performing  the  phasor  summation  at  P.  To  determine 
this  modification  for  the  complex  being  considered,  it  is  first  neces¬ 
sary  to  find  the  component  of  the  linear  velocity  which  is  in  the 
direction  of  P.  This  can  be  accomplished  by  considering  Fig.  A-l,  which 
is  a  projection  of  the  coordinate  geometry  in  the  0  =  90°  plane.  The 
tangential  velocity  of  the  scattering  complex  at  0",  relative  to  the 
primary  origin  0,  is 


109 


•  v 

l 

■ 

b. 

f! 


h 


u 


c 

1"  4 
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The  component  of  this  velocity  which  is  in  the  direction  of  P  is 

vR  =  vT  cos  ? 

=  o)yd  cos  (4>— <J> *  +  ir/2) 

=  -w^d  sin  (<f>— T )  -  (A-3) 

From  the  well-known  Doppler  formula 

A 

a*  .  w  [1  -  i-!-£]  ,  (A-4) 

o  c 

where  c  is  the  speed  of  light  in  free  space  and  r  is  directed  from  the 
point  of  observation  to  the  moving  source,  it  is  seen  that 

;sU0u  +  \  (a-5) 

Rearranging  in  terms  of  wavelength, 

X'  =  X_[l  +  — ]_1.  (A-6) 

°  c 

Equation  (A-6)  given  the  Doppler-shifted  wavelength  X'  in  terms  of  the 
unshifted  value  \Q  and  the  component  of  velocity  of  the  complex  which 
is  towards  P. 

The  relative  unshifted  spatial  retardation  phase  is 

a  =  -kQr  (A— 7) 

where  kQ  ==  2tt/ Xq  and  r  is  the  one-way  path  from  0"  to  P.  Similarly, 
the  shifted  phase  is 


a’  =  -kr. 


(A- 8) 


Thus,  the  apparent  change  in  phase  due  to  the  Doppler  shift  is 


Aa  ■  a  -  a' 

=  r(k’-  k0) 


2^R 

AqC 


(A-9) 


This  result  is  used  in  the  model  simply  by  making  the  substitution 


2irr 


o  .-*■  a  +  sin(d>-<f>‘) 


(A-10) 


for  the  phase  contribution  of  the  complex  being  considered. 

The  proper  coordinate  and  phase  modifications  for  roll  and  pitch 
are  more  complex  because  rotation  about  these  axes  involves  both 
angular  coordinates.  However,  if  a  rotation  of  coordinates  is  first 
performed  to  align  the  z-axis  with  the  roll  or  pitch  axis,  there  is 
once  again  only  an  azimuth  variation.  For  roll,  the  transformations 
are 


<j>"  =  arctan  (-tan  <j>  sin  <j>) 

6"  «  arccos  (sin  0  cos<j>) .  (A-ll) 


After  modifying  the  new  <j>"  for  the  roll  rate  in  a  manner  similar  to 
(A-l) ,  the  new  coordinates  of  P  are  found  from  the  reverse  transforma¬ 
tions 


<j>’  =>  arctan  (tan  0"  sir.  <}>") 

9*  »  arccos  (-sin  9"  cos  (A-12) 


The  analogous  transformations  for  pitch  are 
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FORTRAN  IV  IMPLEMENTATIONS 


This  appendix  presents  flow  diagrams  and  listings  of  the  more 
important  computer  programs  and  subroutines  used  in  the  glint  study. 
Most  of  the  programs  were  originally  prepared  for  an  IBM  7040  computer, 
and  subsequentially  converted  for  use  with  an  IBM  System  360,  Model  50. 

Many  different  TRAJ,  AMPPHI,  and  OUTPUT  subroutines  were  written 
for  use  in  making  various  data  runs;  the  ones  listed  are  illustrative 
of  the  general  programming  techniques  employed.  The  programs  which 
were  used  to  analyze  the  data  and  to  generate  the  various  machine  plots 
are  not  included  in  the  listings. 


Fig.  3-1. — Flowchart  for  main  glint  program,  continued. 
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C  THIS  PROGRAM  COMPUTES  GLINT  ERROR  ALONG  A 
C  PRESCRIBED  TRAJECTORY*  PROVISION  IS  MADE 
C  FOR  MOTION  OF  THE  TARGET. 

C 

IMPLICIT  REAL*8( D) 

COMPLEX*l 6  Z 
REAL  LAMBDA 

REAL*8  PTH<100),P<5),E(5) 

COMMON/COMl/AZI ICO) tEL< 100 > , RD( 100 ) ,NO< 100) , 

1  ROLL  .PITCH, YAW 

COMMON/COM2/AZM,ELV,RG,AMPUOO),PHI  (100)  ,AMPL(100) 
COMfaON/COM4/N, W AVE , CONV  , RM 
COMMON/ COM 5/ AZSI 100) ,ELS( 100) 

C0MM0N/C0M6/ ICNT, JK 
C0MM0N/C0M7/IDC0N, L» LAMBDA tRDI A 
DIMENSION  COMMNTI 20 ) 

C 

PI-3.141593 

C0NV=1B0./PI 

ICNT-0 

C 

C  READ  WAVELENGTH,  TIME  INCREMENT,  MINIMUM 
C  RANGE,  APERTURE  DIAMETER,  ROLL  PITCH  AND  YAW 

C  RATES,  ELV  ERROR  SUPPRESS  CODE,  COORDINATES  OF 

C  SCATTERING  COMPLEXES,  AND  COMMENT  CARD 

READ(5,100)LAMBrA,TINC,RMIN,RDIA,ROLL, 

1  PITCH, YAW, NOELV 
READ{5,101) IDCON,N 
DO  5  NN-1 » N 

READ! 5, 102 ) I , RDI I ) , AZ< I ) ,EL( I ) 

AZ( I )=AZ( I )/CONV 
EL( I  )=EL( I )/CONV 
5  NO(NN)»I 

READ<5,103)(C0MMNT(KK),KK=1»20) 

C 

WAVE=(2.*PI)/LAMBDA 
DPOI A=WAVE*RDI A 
D2Pl-2.*PI 
DCONV-CONV 

CALL  TRAJ(L,0.,RG, ELV, AZM, CONV) 

ELV=ELV*CONV 
AZM=AZM*CONV 
WRITE (6, 200) IDCON,L 
WR I T  E ( 6  »  COMMNT ) 

CALL  AMPPHI 
WR ITE (6, 206) JK, RM 
WRITE(6,201)RG,ELV,AZM 
WRITE (6, 202) ROLL, PITCH, YAW 


Fig.  B-2. — Listing  of  main  glint  program. 


ooo  oo  oo  o  o  ooo  or*  ooo  oo 
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WRITE(6,203) LAMBDA,RD1A 
C 

C  INCREMENT  TIME 

DO  90  K=1 ,  1001 
T=(K-U*TINC 

CALL  FIELO  POINT  COORDINATES  FOR  THIS  TIME 
GALL  TRAJ( L*T,RG»ELV»AZM»CGNV) 
IF(RG.LT.RMIN)  GO  TO  97 
ORG-RG 


MODIFY  COORDINATES  TO  INCLUDE  TARGET  MOTION#  IF 
REQUIRED 

IF<  ( ROLL+P ITCH* YAW ) .GT.O. I CALL  MOTION! T,AZM, 

1  ELY* ROLL, PITCH, YAW) 

BEGIN  FINDING  AMP  AND  PHASE  OATA  FOR  COMPLEXES 
00  15  NN=1,N 
I=N0(NN) 

CALL  DI ST ( AZM, ELV, RG, AZ<  I ) ,EL ( I) »RD( I ) , PTHI I ) ) 
RDIST=PTH( I ) 

COMPUTE  AZS  AND  ELS— THE  EQUIVALENT  AZM  AND  ELV 
IN  THE  ITH  COORDINATE  SYSTEM 
X=RG*SIN( ELV ) *C0S(  AZM ) 

1  -RDU)*SIN<FLU>)*COS(AZ(I)) 

Y=RG»COS( ELV )-RD( I ) »C0S(EL( I ) ) 

ELS  I  I )* ARGOS (Y/RDI ST ) 

AZS(  I)  =  ARC0S{X/(RDIST*SIN(ELSCI ))) ) 

CALL  AMP  AND  PHASE  VALUES  FOR  ITH  COMPLEX 
15  CALL  AMPPHK  I,  AZ${  I )  ,ELSU  ) ,  AMP  ( I )  ,PHI  (!) ) 

COMPUTE  TOTAL  AMP  AND  PHASE  AT  THESE  COORDINATES 
CALL  SIGNAL! 1,K,DX,DY,PTH) 

Z=DCMPLX { DX, DY ) 

A*CDABS{ Z ) 

DG=DSIGN( DARCOS ( DX/CDABS(Z)),DY) 

COMPUTE  HALF-ANGLFS  SUBTENDED  BY  APERTURE 
RAD=RDI A/2. 

EL3*ATAN( RAD/RG) 

AZ3=ATAN(RAD/(RG*SIN(ELV) ) ) 

BEGIN  COMPUTATION  OF  ERROR  IN  THE  TWO  TRACKING 
PLANES 

DO  29  J  =  1 » 5 

GO  TO  (21, 22,23, 27, 30), J 


Fig.  B-2. --Listing  of  main  glint  program,  continued. 
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21  AZM=AZM+AZ3 
60  TO  28 

22  AZM* AZM-2 • *AZ3 
GO  TO  28 

23  DC*P( 1> 

0F=P! 2) 

CORRECT  FOR  ANY  PHASE  CYCLE  CHANGE  ACROSS  APERTURE 

24  IF  ! DABS I DC-DG ) .IT . . 50D0. AND. OABS I DG-DF ) • LT. . 50D0 ) 
1  GO  TO  26 

I F I (DC-DG)*! DG-DF) .GE.O. )  GO  TO  26 
IFtDC.GT.DF)  GO  TO  25 
DC*DC+D2PI 
GO  TO  26 

25  0C=DC-D2PI 

COMPUTE  SINGLE-PLANE  ERROR 

26  E!J)*!DRG*!DC-DF))/DPDIA 
IF( J.EQ.5 )  GO  TO  32 
IFINOELV.EQ.O)  GO  TO  31 
AZM*AZM+AZ3 
ELV-ELV+EL3 
GO  TO  28 

27  ELV»ELV-2.*EL3 

28  CALL  SIGNAL!  2,  Ki,DXf.OY,PTH) 

Z*OCMPLX(DX«DY) 

29  P ! J ) =DS I GN  < DARCOS ( DX/CDABS ! Z ) ) »  OY ) 

30  DC«P!3) 

DF=P!4) 

GO  TO  24 

31  E!5)*0. 

COMPUTE  TOTAL  ERROR 

32  ERRQR«DCONV*OAT  AN2 ! DSQRT! E ! 3 )**2+E ( 5 ) **2 ) , DRG) 
G*DG 

C 

C  TRANSFER  RESULTS  TO  OUTPUT  OEVICES(S) 

90  CALL  OUTPUT! K» T ,E(3),E!5) , ERROR , A,G, AMPL ) 

97  CONTINUE 

C 

100  FORMAT! ?F10.0»I1) 

101  FORMAT! 212) 

102  FORMAT! 12, 3F10.0) 

103  FORMAT! 20 A4) 

200  FORMAT! 1H1,34HGLINT  ERROR  FOR  CONFIGURATION  NO.  , 
1  I2»/»22H  ALONG  TRAJECTORY  NO.  ,13) 

201  FORMAT! 1H0,29HINITIAL  SEEKER  POSITION— RG=  , 

1  F7. 1, 2H  M,/f 26X,4HEL*  ,F7.3,4H  DEG,/, 


Fig.  B-2. — Listing  of  main  glint  program,  continued. 
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2  26X»4HAZ=  ,F7.3,4H  DEG) 

202  FORMAT ( 1H0» 21HT ARGET  MOTION — ROLL=  ,F8.3, 

1  8H  RAD/SEC, /,16X,7HPITCH=  ,F7.3,8H  RAD/SEC,/, 

2  16X , 5HYAW=  ,F9.3,8H  RAD/SEC) 

203  FORMAT J 1H0* 12HWAVELENGTH=  ,2PF11.3^3H  CM,/, 

1  13H  SEEKER  OIA=  ,F11.3,3H  CM//) 

206  FORMAT ( 1H0 , *  AMP PH I  DATA  SOURCE  NO.  ',13, 

1  •  AT  RM  =  * ,F7. 1* •  M') 

STOP 

END 


Fig.  B-2. — Listing  of  main  glint  program,  continued. 
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SUBROUTINE  SIGNAL<L,k»DX*DY,PTH) 

C 

C  THIS  SUBROUTINE  COMPUTES  THE  TOTAL  RECEIVED 
C  SIGNAL  FROM  A  GROUP  OF  SCATTERING  COMPLEXES. 

C  PROVISION  IS  MADE  FOR  DOPPLER  DUE  TO  MOTION 
C  OF  COMPLEXES. 

C 

IMPLICIT  REAL*8  (D) 

REAL*8PTH{ 100) 

C0MM0N/C0M1/AZ<  100) ,  EL  <100 ) ,  RD(  100 )  ,  NOUOO) , 

1  ROLL, PITCH, YAW 

C0MM0N/C0M2/AZM, ELV»RG» AMP< 100) , PHI ( 100) ,AMPL( 100) 
COMMON/COM3/DCONV, DWAVE ,DPTH,DPHI, I 
C0MM0N/C0M4/N, WAVEfCONV.RM 
C 

IFUL+K)  .GT.2)  GO  TO  70 
DREDU=3.6E7/(WAVE*C0NV) 

DPMAX-1. 3D0*DREDU 
OCONV=CONV 
DWAVE=WAVE 
DRM=RM 
C 

C  BEGIN  COMPUTATION  OF  TOTAL  RECEIVED  SIGNAL 
70  DX=0. 

DY=0. 

DO  80  NN*=  1 ,  N 
IsNO(NN) 

C 

C  FIND  PATH  OF  PROPAGATION  ASSOCIATED  WITH  COMPLEX 
IF  <  L • EQ. 1 )  GO  TO  72 

CALL  DIST(AZM,ELV,RG,AZU),EL<  I  )  *RD(  I  )  ,DPTH ) 

GO  TO  73 

72  DPTH=PTH(I) 

C 

C  MODIFY  PHASE  TO  INCLUDE  DOPPLER,  IF  REQUIRED 

73  DPHI  =  PHI  ( I  ) 

I F < ( ROLL+P I TCH+YAW ) .GT.O.  )  CALL  OOP 
DPATH=PTH( I )+DPTH-2.D0*DRM 
PATH1=DPATH 
C 

C  REDUCE  PATH  TO  LESS  THAN  130,000  LAMBDA 

77  IF ( DPATH.LT .DPMAX )  GO  TO  78 
DPATH=DPATH-OREDU 

GO  TO  77 
C 

C  FIND  RECEIVED  SIGNAL 

78  E=AMP ( I ) / 2  0 • 

DAMPL=4 10. **E ) /PATH1 


r 

IU 

Fig.  B-4. — Listing  of  subroutine  SIGNAL, 
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AMPL(NN)=DAMPL 

DPHSE=<DPHI/DCONV)-(OWAV£»DPATH) 
OX1=OAMPL*DCOS( DPHSE ) 
OY1=OAHPI*DS IN<  DPHSE ) 

DX=DX+DX1 
60  DY=DY+DY1 
RETURN 
ENO 


Fig.  B-4.— Listing  of  subroutine  SIGNAL,  continued. 


COMPUTE  POLY 
ARGUMENT,  X 


CALL  POLY 


— Flowchart  of  subroutine  BIST 


o  o  o  o  o  o 


SUBROUTINE  DIST ( AZM.ELV , RG. AZ I t£LI i RDI tDPATH) 


THIS  SUBROUTINE  COMPUTES  THE  1-WAY  PATH  BE1WEEN 
THE  FIELD  POINT  WITH  SPHERICAL  COORDINATES 
(RGt ELV» AZM)  AND  THE  POINT  WITH  COORDINATES 
I RDI »ELI t AZ I ). 

IMPLICIT  REAL*8  ( D ) 

DAZ=AZM 
DEL=ELV 
DRG=RG 
DA I  =  AZ  I 
DEI =£L I 
DO  I =RD I 
C 

DX=DCOS( DEL) *DCO$ ( DE I )+DS I N( DEL ) *OS l N(DE I ) 

1  *DCOS(  DAZ.-DAI ) 

C 

DS=I .DO/DRG 
DO  50  K=1 f 9 
DR=D0I**K/DRG**{K+1) 

IF (DR.LT. 1 .0-2 C )  GO  TO  10 
CALL  POLY ( K » DX  t  DP ) 

OS=OS+DR*DP 
50  CONTINUE 

WR I TE { 6  »600 ) 

600  FQRMATdHO,  ’SERIES  NOT  CONVERGING') 

10  DP A  TH= 1 . DO/DS 

RETURN 
END 


Fig.  B-6. — Listing  of  subroutine  DIST. 


o  o  n  o 


. . .SUBROUTINE  PULY l L » DX ,QP ) . . . . . . .  . 

. THIS.. SUB  ROUTINE  COMPUTES  THE  LEGENDRE  POLYNOMIAL 

OF  ORDER  L  AND  ARGUMENT  DX. 

IMPLICIT  REAL*8  (D) 

.  GO  TO  ( I  *  2  » 3.*.4.f  5  » 6 » 7 ,8 . 9 )  .».L . 

1  DP=DX 
. GO  TO  10 

2  DP  =  b.D-14(3.D0*DX5!t*2-l.D0) 

GO  TO  10 

3  DP  =  5.D-l’M5.D0*DX*#3-3.D0*tDX) 

. GO  TO  .10. 

4  DP  =  1.2  5D-1=M3.5D1*DX**4-3.D1*DX**2+3.D0) 

. GO  TO.  10  .  . . . - . . 

5  DP  =  1 .2  5D- 1* ( 6. 3D1*DX**5-7.D1*DX**3+1.5D1*DX) 

GO  TO  IQ  . - . - . . 

6  DP =6 .250-2* (2.3ID2*DX**6-3. 15D2*DX**4 

. I  +1*05D2*0X**2-5.D0) . . 

GO  TO  10 

7  DP=6.25D-2*(4.29D2*DX**7-6..9.3D2*0.X**5. . . 

1  >3. 1502*DX**3“3.5D1*DX) 

GO  TO  10 

8  DP=7.8125D-3*(  6. 43 503*DX**8- 1 . 20  1 2D4*DX**6 

1  +6. 93D3*DX**4-1.26D3*DX**2+3. 5D11 . 

GO  TO  10 

9  DP=7.8125D-3*(1.2155D4*DX**9-2.574D4*DX**7 

1  + 1 . 801 8D4*DX**  5-4 .6203*0 X**  3+3 . 15D2*DX ) 

10  CONTINUE  . 

RETURN 
END 


Fig,  B-7.—  Listing  of  subroutine  POLY 
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SUBROUTINE  MOTION! T,AZM,ELV, ROLL, P ITCH, YAW ) 

C 

C  THIS  SUBROUTINE  CHANGES  THE  ASPECT  ANGLE  TO 
C  ACCOUNT  FOR  TARGET  MOTION.  IT  IS  VALID  ONLY  FOR 
C  ROTATION  ABOUT  THF  Z-AXIS  (YAW)  AND  IS  RESTRICTED 
C  TO  ASPECT  ANGLES  TN  THE  ELV  *  90DEG  PLANE. 

C 

IF( (ROLL+PITCH) .NE.O.) WRITE! 6,200) 

IF! ELV.GT.I.575.0R.ELV.LT. 1.565) WRITE (6, 200) 

200  FORMAT < 1H0, 28HMOT ION  CALLED  WHEN  NOT  VALID) 

C 

AZM=AZM+YAW*T 
DO  5  J=1 , 20 

IFIAZM.LT. 6. 2831853)  GO  TO  6 
AZM=AZM-6. 2831853 
CONTINUE 
CONTINUE 
RETURN 
END 


Fig.  B~9. — Listing  of  subroutine  MOTION. 


DEFINE 

CONSTANTS 


Fig.  B-10. — Flowchart  of  subroutine  DOP. 
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SUBROUTINE  OOP 
C 

C  THIS  SUBROUTINE  INCLUDES  DOPPLER  SHIFTS  INTO  THE 

C  PHASE  FROM  A  GIVEN  COMPLEX*.  IT  IS  VALID  ONLY 

C  FOR  A  ROTATION  ABOUT  THE  Z-AXIS  (YAW)  AND  IS 

C  RESTRICTED  TO  ASPECT  ANGLES  IN  THE  ELV  *  90DEG 
C  PLANE. 

C 

IMPLICIT  REAL*8(D) 

COMMON/COM 1/ AZ (100 ) , EL  <100 ) , RD ( 100 ) , N0( 1 00 ) , 

1  ROLL, PITCH, YAW 

COMMON/ COM3/DCONV ,  DWAVIE « DPATH, OPHI  ,1 
C0HM0N/C0M5/AZS< iOO) ,ELS( 100) 

C 

IF< ( ROLL+P ITCH ).NE.O.) WRITE!  6, 200) 
IF(ELS(I).GT*1»575.0R«ELS(I).LT. I. 565) WRI TE (6,200) 
200  FORMAT! 1H0,*DQP  CALLED  WHERE  NOT  VALID*) 

C 

OYAW*YAW 
DIaRD( I ) 

DAZ I=AZ ( I ) 

DAZS*AZS( I ) 

DPHI*DPHI-{ ( DPATH*DYAW*DCONV*DWAVE*DI ' 

1  *DSIN< DAZI— DAZS) )/3.D8) 

RETURN 

END 


Fig,  B-ll. — Listing  of  subroutine  DOP. 
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SUBROUTINE  OUTPUT (K»T»EAZ* EEL* ER» A* P»AHPL) 

THIS  SUBROUTINE  WRITES  THE  PROGRAM  RESULTS  INTO 
MAGNETIC  DISK  STORAGE*. 

REAL  LAMBDA 

COMMON/COMl/DUMt 400), ROLL, PITCH, YAW 
COMMON/COM4/N, WAVE , CONV,RM 
C0MM0N/C0M6/ ICNT, JK 
COMMON/COM7/ IDCON, L  »  LAMBDA, D I A 
DIMENSION  AMPL(IOO) 

IFIK.GT.1)  GO  TO  10 

WRITE  PARAMETERS 
NUMBER*1 

WR I TE ( 8  > NUMBER*  TDCON ,L»JK, KM , LAMBDA , D I A , ROLL  * 
1  PITCH, YAW 
WRITE!6*205) 

WRITE  A  LINE  OF  RESULTS 
10  KRITE(8)T*  EAZ, EEL, P, A# ( AMPL ( J) »  J=1 ,4) 

WRITE! 6, 204 IT, EAZ»EEL»ER,A,P* ! AMPL (J) » J*i»4) 
IF! K.EQ. 1001 )  END  FILE  8 

204  FORMAT ! 1H  ,F9.3,2F10.4*F9.3, 7E13.4) 

205  FORMAT! 1H  ,4X,4HT1ME,4X,3HEAZ,7X,3HEEL,.3X, 

1  ‘TOT  ANG  ERRS5X,  *AMP«t.lOX, 'PHASE* ) 

RETURN 

END 


Fig,  B-12. — Listing  of  typical  subroutine  OUTPUT. 
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Stja.iuUTl.Nifc-  7^J(L-*T*RG*CLV^A4M*£CNV)  -  - 

TRAJECTORY  NO.  11  -------  -  - 

THIS  SUBROUTINE  GENERATES  A  SET  Ofc  H  fci.0 

CGOkO I  NAT  tS .  AS  A  FUNCTION  Or  TIMt,  WHICH  OESCRU3E 

A  Dfcf.IkCO  TKAUfcCTORY  . . .  -  -  - 

L=  1 1  ------  ------- 

RG=4650,-3.7  5*T 

4 Z m=4S. /G ONV+A T^AN2 (I  50.  *Sl  Nl 6.23 3 2 *T RG-SGCU 4 / 

1  3760. }»RG) 

EL  V=9().  /C0NV+ATAN2  ( 5C.0*S I  NC&-»-2-6'32-*4 RG--9CG.  >/ 

1  3  750. ) »  RG l 

RETURN . -  -  -  --  . - . 

tNO 


b 


Fig.  B-13. — Listing  of  typical  subroutine  TRAJ 


SUBROUTINE  AMOBMt  (  I  ,AZS,ELS,AMP, PHIl 


C 

C  A“PPWI  DATA  SOURCE  NO.  QQ 
C 

C  THIS  SUBROUTINE  SYNTHESIZES  AMP  AND  PHASE  VALUES 
C  FOR  COMPLEXES  20.  21.  22,  AND  23 
C 

C0MM0N/C0M4/N,  WAVE ,CONV*RM 
C0MM0N/C0M6/ICNT » JK 
IF  I ICNT.NE.O)  GO  TO  P9 
JK=99 
R  M  =  0  « 

I CNT  =  1 

PI = IRQ. /CONV 
GO  TO  PS 

89  CALL  RDUND( AZS.E  LS.AZ.EL) 

K= I-  19 

GO  TO  (5,10,15,15) , K 

5  A  =  (SIN(FU*SIN<AZ)*COS(380.*AZ)I**2  . —  - 

PHI=180.*C0S(720.*AZ) 

GO  TO  ?0 

10  A=( .9*SIN(EL)*SIN(AZ)*COS(720.*AZ) )**2 

DHI=180.*C0$(360.*AZ)  . 

GO  TP  20 

15  IF(AHS(  AZ)  .LE.PI/2.)A=(SIMEU*C0S(AZ) 

1  *COS( 360.*AZ))**2 

IF(A8S(AZ).GT.PI/2.)  A=C. 

PHI=130.*C0S(7?0.*AZ) 

20  IF(A.EO.O.)  A=1 . E-50 

AMP=10.*ALOG10(A) 

88  CONTINUE 

RETURN 
END 


Fig.  B-14 — Listing  of  typical  subroutine  AMPPHI. 


SUBROUTINE  ROUNDC A,E,AZ,EL ) 


C 

C  THIS  SUBROUTINE  ROUNDS,  TO  THE  NEAREST  'R',  THE 

C  VARIABLES  A  ANO  E.  THE  RESULTS  ARE  AZ  AND  EL. 

C 

C  THE  DESIRED  ROUNDTNG  FACTOR,  IN  RADIANS,  IS 

C  PUNCHED  INTO  THE  NEXT  CARD 
R=. 01745 

I F ( R.EQ.O. )  GO  TO  90 
E=-l.— AINT( ALOGIO(R) ) 

FACT=10.**E 
RND=R*FACT 
AZ=ABS( A J*FACT 
AZ=RN0*AINT ( (AZ/RNDI+.5) 

AZ=SIGN( (AZ/FACT),A) 

EL=ABS( E ) *FACT 
EL=RND*AINT( (EL/RNDJ+.5) 

EL*SIGN( (EL/FACT), E) 

GO  TO  91 

90  AZ= A 
EL=E 

91  CONTINUE 
RETURN 
END 


Fig.  B-15. — Listing  of  subroutine  ROUND. 
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